AD-A069  291  NAVAL  ACADEMY  ANNAPOLIS  MD  DIV  OF  ENGINEERING  AND  WEAPONS  F/G 
SYNTHESIS  OF  OPTIMAL  LADDER  NETWORKS. (U) 

NOV  78  T S LIM 

UNCLASSIFIED  USNA-EW-16-78  NL 


F/G  12/2 


1 OF  2 

A064  29! 

■ 

IB. 

• 

■ ' 

- 

* 

- 

3 

3 

- 

UNITED  STATES  NAVAL  ACADEMY 
DIVISION  OF 

ENGINEERING  AND  WEAPONS 

ANNAPOLIS,  MARYLAND 

D D 

~ DISTRIBUTION  STAT . " 

Approved  for  public  ; . 

Distribution  Unliro: 


idFT"  u 


DEPARTMENT  OF  THE  NAVY 

UNITED  8TATE8  NAVAL  ACADEMY 
ANNAPOLI8,  MARYLAND  21402 


DIVISION  OF  ENGINEERING  AND  WEAPONS 


Report  EW-16-78 

Synthesis  of  Optimal 
Ladder  Networks 


m 

•X 


Tian  S.  Llm* 
November  1978 


1 tttailM  Mr  | 

V* 

mnt 

•M 

■•n  mum  'n 

M»JKMM6S9 

JtflTIffttTI 

a 

n 

liSTIIBSTIOH, 

cost* 

Bt».  Z1 

: nt  v :*tcuu. 

Plj 

! 


♦Assistant  Professor 
Electrical  Engineering  Department 
U.  S.  Naval  Academy 
Annapolis,  Maryland 


D 


i 


D 

FE3  5 <ST9 

JljElLL.. 


it  i 

» l 
» j 

• i 


This  report  was  prepared  in  connection  with  research  grants  from 
the  U.  S.  Naval  Academy  Research  Council  (NARC). 


DISTRIBUTION  STATES--:?-:  \ \ 

Approved  for  public  ici'x-. 
Distributicn  UnJ 5 milt'd 


ItMHIUT  (.LHHintHTHW  or  THU  r<8l  PWIM  BjM  «W«m 


omnnun  i.i 


REPORT  DOCUMENTATION  PAGE 


NUMB 

EW-16-78 


4.  TITLE  fond  Subtlllm  1 


NTHESIS  OF  OPTIMAL 
DOER  NETWORKS  , 


LADDER  NETW0R1 


AUTHOR/.J 


mm 


United  States  Naval  Academy 
Annapolis,  Maryland  21402 


II.  CONTROLLING  office  name  and  adoress 

United  States  Naval  Academy 


Hll  lil  H LI*/'1 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT'S  CATALOG  NUMBER 


14.  DISTRIBUTION  STATEMENT  (at  (hi*  R*par() 


NAME  k ADDRESS (II  different  from  Controlling  Office) 

IS.  SECURITY  CLASS,  (ol  thio  report) 

Unclassified 

| 

:CL  ASSI FIC  ATI  ON/  OOWN  OR  AOIN  0 

hedule 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  obetroet  entered  In  Block  20,  If  different  from  Report) 


IS.  KEY  BOROS  (Continue  on  reretee  old*  It  mcnaor  and  Identity  by  4 lock  number) 


Ladder  networks,  optimal 


tO.  ABSTRACT  (Continue  on  roeereo  o!4o  If  neeeeoory  end  Identify  by  block  number) 

This  project  treats  synthesis  procedures  for  optimal  ladder  networks. 
The  first  part  of  the  project  deals  with  a chain-matrix  decomposition 
technique  for  realizing  two-port  LC  ladder  networks.  Two  chain-matrix 
decomposition  theorems  are  proved.  These  theorems  state  the  necessary 
and  sufficient  conditions  that  must  be  satisfied  by  a ^2x2 J matrix  Sn(s) 
that  is  to  be  decomposed  into  a product  of  simple  matrices,  i.e.,  ^ ^ ■ 


S„(s)  - K,K2. 


It  is  found  that  the  zeroes  of  the  adjacent  elements 


I JAN  T» 


EOinOM  OF  I NOV  BB  It  OBSOLETE 
B/N  0^}02Y0I4*  4601  I 


SECURITY  CLASSIFICATION  OF  THIS  FAOC  ("bon 


A b 


' b -v.  i 


-UCUHITY  CLASSIFICATION  OF  THI»  PAOCOW.«n  Pl«  lnl»l»<)  / 

of  Sn(s)  aUemate  pairwise  along  the  jw  axis  in  the  s-plane 


It  is 


also  found  that  if  the  matrix  Sn(s)  is  the  overall  chain  matrix  of  an  LC 
ladder  network,  then  each  of  the  simple  matrices  Ki(s)  for  1 3 l,2,*’*n 
represents  a simple  LC  ladder  section.  (_ 


• *> 


k 


These  decomposition  techniques  are  then  applied  to  the  design 
of  filters.  Of  particular  interest  are  Butterworth,  Chebyshev,  and  Bessel 
filters  with  single  and  double  terminations.  These  filters  are  designed 
by  the  decomposition  of  chain  matrices  whose  elements  are  predetermined 
by  the  orthogonal  polynomials  that  approximate  the  Ideal  filter 
characteristics.  , 

The  synthesis  of  inhomogeneous  ladder  networks  is  also  Investigated. 
By  means  of  the  mapping  p = y(s)z(s),  the  optimal  process  for  RC  ladders 
Is  generalized  to  synthesize  optimal  Inhomogeneous  ladders. 

The  last  part  of  this  project  treats  the  synthesis  of  double- 
terminated  optimal  RC  ladders  whose  eigenvalues  follow  certain  patterns. 
Several  formulas  are  found  to  facilitate  the  calculation  of  element 
values  for  optimal  ladders  whose  eigenvalues  are  pj  3 ip]  as  well  as  for 
those  ladders  whose  eigenvalues  are  pi  = 12pi.  These  formulas  reduce 
the  tedious  process  of  optimal  synthesis  to  simple  formula  substitutions. 


SSCUMTY  CLASSIFICATION  OF  THIS  FAOSnVMn  DM* 


r 


IV 


ABSTRACT 

This  project  treats  synthesis  procedures  for  optimal 
ladder  networks.  The  first  part  of  the  project  deals  with  a 
chain-matrix  decomposition  technique  for  realizing  two-port 
LC  ladder  networks.  Two  chain-matrix  decomposition  theorems 
are  proved.  These  theorems  state  the  necessary  and  suffi- 
cient conditions  that  must  be  satisfied  by  a [2x2]  matrix 
Sn(s)  that  is  to  be  decomposed  into  a product  of  simple 
matrices,  i.e.,  S fs)  * K, K,***K  . It  is  found  that  the 
zeroes  of  the  adjacent  elements  of  Sn(s)  alternate  pairwise 
along  the  jw  axis  in  the  s-plane.  It  is  also  found  that  if 
the  matrix  S fs)  is  the  overall  chain  matrix  of  an  LC  ladder 
network,  then  each  of  the  simple  matrices  K^(s)  for 


i = 1,2, •••n  represents  a simple  LC  ladder  section. 

These  decomposition  techniques  are  then  applied  to  the 
design  of  filters.  Of  particular  interest  are  Butterworth, 
Chebyshev,  and  Bessel  filters  with  single  and  double  ter- 
minations. These  filters  are  designed  by  the  decomposition 
of  chain  matrices  whose  elements  are  predetermined  by  the 
orthogonal  polynomials  that  approximate  the  ideal  filter 
characteristics . 

The  synthesis  of  inhomogeneous  ladder  networks  is  also 
investigated.  By  means  of  the  mapping  p = y(s)z(s),  the 
optimal  process  for  RC  ladders  is  generalized  to  synthesize 
optimal  inhomogeneous  ladders. 
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The  last  part  of  this  project  treats  the  synthesis  of 

double- terminated  optimal  RC  ladders  whose  eigenvalues 

follow  certain  patterns.  Several  formulas  are  found  to 

facilitate  the  calculation  of  element  values  for  optimal 

ladders  whose  eigenvalues  are  p^  = ipj  as  well  as  for  those 

2 

ladders  whose  eigenvalues  are  p^  = i pj.  These  formulas 
reduce  the  tedious  process  of  optimal  synthesis  to  simple 
formula  substitutions. 
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CHAPTER  1 


REVIEW  OF  THE  LITERATURE  AND  SCOPE  OF  THE  THESIS 

The  study  of  ladder  network  synthesis  began  in  the 
1920s.  The  first  paper  dealing  explicitly  with  the  realiza- 
tion of  a one-port  whose  impedance  is  a prescribed  function 
of  frequency  was  written  by  Cauer  [6]  in  1926.  Cauer's 
method  is  based  on  continuous  fraction  expansions.  Cauer 
also  introduced  one  of  the  most  direct  methods  for  synthesiz- 
ing two-port  networks  resulting  in  ladder  structures  [7]. 
Since  then,  a number  of  synthesis  procedures  for  two-port 
networks  have  been  discovered  [24].  However,  very  few  of 
these  methods  made  use  of  chain  parameters. 

One  of  the  earlier  works  on  the  synthesis  of  two- 
elemcnt-kind  two-port  networks  using  chain  matrix  factoriza- 
tion was  published  by  Hunt  [10].  In  Hunt's  method,  a given 
chain  matrix  can  be  factorized  as  a product  of  elementary 
matrices,  each  representing  one  branch  of  a section  of  a 
ladder.  However,  Hunt  did  not  give  any  realizability  con- 
dition on  the  overall  chain  matrix.  Therefore,  at  the  end 
of  each  step  of  Hunt's  procedure,  one  cannot  be  sure  of  the 
realizability  of  the  remaining  matrix.  The  procedure  tends 
to  be  one  of  trial  and  error.  Lee  and  Brown  [15]  were  the 
first  to  prove  the  necessary  and  sufficient  conditions  of 
the  overall  chain  parameters  for  the  realization  of  two- 
element-kind  ladder  networks  using  chain  matrix  decomposition 
techniques . 


* 
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The  motivation  of  this  thesis  is  due  in  part  to  the 
resurgence  of  interest  in  the  optimal  realization  of  two- 
el  ement  - k i nd  two-port  networks  with  specified  transmission 
zeros.  Optimal  synthesis  of  RC  ladder  networks  with  given 
terminations  was  first  introduced  by  Kuh  [12]  in  1958.  He 
realized  a grounded  RC  two-port  from  a given  transfer 
function  with  distinct  poles  on  the  negative  real  axis, 
lie  used  the  Lagrange  multipliers  method  to  obtain  the 
optimization  and  the  resulting  optimal  network  has  maximum 
gain.  integrated  circuit  technology  spurred  a revival  of 
interest  in  optima]  synthesis  of  RC  ladders  nearly  two 
decades  after  Kuh’s  discovery.  Optimal  synthesis  of  RC 
ladder  networks  using  chain  parameters  in  the  s-domain  was 
discovered  by  Protonotar ios  [21]  in  1974.  He  formulated 
the  realization  technique  of  RC  ladders  with  n shunt 
capacitors  and  n+1  series  resistors.  His  realization  is 
obtained  from  a given  chain  parameter,  B(s)  , resulting  in 
a network  having  minimum  total  shunt  capacitance.  Pro- 
tonotarios  also  formulated  the  synthesis  of  RC  ladders  with 
n shunt  capacitors  and  n series  resistors.  This  realization 
is  derived  from  a given  chain  parameter,  A(s) , resulting  in 
a network  having  minimum  total  resistance  - total  capacitance 
product.  Almost  simultaneously,  Stein  and  Salama  [23] 
developed  optimal  synthesis  procedures  for  the  same  RC 
ladders  and,  furthermore,  showed  that  the  optimal  RC  ladder 
with  n+l  resistors  is  symmetrical  and  the  optimal  RC  ladder 
with  n resistors  is  antimetrical . 
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Chapter  2 of  this  thesis  deals  with  lossless  two-port 
synthesis  of  ladder  networks,  using  the  chain  matrix  de- 
composition method  which  is  an  extension  of  the  work  done 
by  Lee  and  Brown.  Two  chain  matrix  decomposition  theorems 
are  proved.  These  theorems  give  the  necessary  and  suffi- 
cient conditions  of  the  chain  parameters  for  the  existence 
of  cascaded  lossless  ladder  networks.  The  results  of  these 
decomposition  theorems  and  synthesis  algorithms  are  applied 
to  the  design  of  filters  in  Chapter  3.  The  first  half  of 
Chapter  3 is  concerned  with  filter  networks  that  are  ter- 
minated in  a one-ohm  load  resistor  [17],  The  second  half 
of  Chapter  3 deals  with  double-terminated  filter  design 
using  chain  parameter  technique  [18].  Classical  filters 
such  as  Butterworth,  Chebyshev,  and  Bessel  are  synthesized 
through  the  use  of  the  decomposition  of  chain  matrices. 

It  was  found  that  with  suitable  terminations  these  filters 
are  optimal. 

Chapters  4 and  5 extend  the  work  done  by  Lee  and  Brown 
on  inhomogeneous  ladder  networks  [16]  as  well  as  the  work 
done  by  Kuh,  Protonotarios , and  Stein  on  optimal  RC  ladders. 
Chapter  4 deals  with  optimal  synthesis  of  the  inhomogeneous 
ladder  networks  as  shown  in  Figure  4.1  and  Figure  4.2.  By 
means  of  the  mapping  of  p * y(s)z(s),  the  optimal  process 
was  generalized  to  synthesize  inhomogeneous  ladders. 

Chapter  5 is  concerned  with  synthesis  procedures  of  certain 
classes  of  optimal  RC  ladder  networks  as  shown  in  Figure  5.1 
with  Rj  = Rn+1*  The  state  equation  is  [V]  = [A] [V]+[B]u 
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where  the  state  variable  [V]  represents  the  voltages  across 
the  capacitors,  u is  the  forcing  function,  and  [A]  and  [B] 
are  constant  matrices  of  a network  of  given  topology  and 
components.  It  was  found  that  if  the  tridiagonal  matrix 
f A J has  eigenvalues  pi  « ip^,  then  its  diagonal  elements 

are  equal  and  each  is  equal  to  the  sum  of  eigenvalues 
divided  by  the  number  of  eigenvalues.  Several  formulas 
have  been  found  and  their  applications  demonstrated.  These 
formulas  reduce  the  tedious  process  of  optimal  synthesis  to 


simple  substitution  in  formulas.  Impulse- response  and  step- 

I resnonse  curves  as  well  as  tahles  r»f  resirl>iec  n-ro  nitron 


CHAPTER  2 


SYNTHESIS  OF  LADDER  NETWORKS  BY  CHAIN  MATRIX  DECOMPOSITION 

Section  2,1  Introduction 

The  object  of  this  chapter  is  to  derive  and  demonstrate 
the  synthesis  procedures  of  the  two-element-kind  lossless 
ladder  networks  by  chain  matrix  decomposition  method. 

Cascade  synthesis  procedures  by  Darlington  and  Youla  are 
found  useful  for  immittance  function  realization  but  their 
methods  are  different  approaches  than  the  chain  matrix 
techniques . 

The  properties  of  entries  of  a two-by-two  matrix  to  be 
synthesized  are  derived  as  lemmas.  Theorems  of  decomposi- 
tion of  this  matrix  into  a product  of  simple  matrices  are 
then  followed.  When  the  ladder  network  is  shown  as  in 
Figure  2.6,  the  zero  distribution  of  each  overall  chain 
matrix  entry  has  a definite  pattern.  Therefore,  suppose 
there  comes  a two-by-two  matrix  with  each  entry  having  this 
same  zero  distribution  pattern,  then  this  matrix  can  be 
used  as  a chain  matrix  of  a realizable,  two-element-kind 
lossless  ladder  network. 

The  necessary  and  sufficient  conditions  of  a matrix 
Sn(s)  which  is  to  be  decomposed  into  a product  of  simple 
matrices  are  obtained  as  shown  in  the  theorems,  and  these 
simple  matrices  K^(s)  for  i * 1,  2,  ...,  N are  the  chain 
matrices  each  of  which  represents  a simple  ladder  network. 

A simple  decomposition  algorithms  is  developed  in 


A 


Section  2.4.  Using  this  algorithms  implemented  by  digital 
computer  programs,  synthesis  of  two -element -kind  lossless 
ladder  networks  can  be  accomplished  with  ease  and  with 
little  expenditure  of  time. 

Section  2.2  Distribution  of  Zeros  of  Chain  Parameters  of 
Cascaded  Ladders 

The  following  five  lemmas  provide  the  properties  of  the 
A,  B,  C,  D entries  of  the  overall  chain  matrix.  These 
lemmas  will  be  used  for  the  proof  of  the  decomposition 
theorems  in  the  next  section. 

rA(s)  B (s) 


Lemma  1 Let  a matrix  S fs)  = 

n 


C (s)  D(s) 


If 


(1)  Det  Sn(s)  = 1, 

(2)  A(s),  B(s),  C(s)  and  D(s)  are  polynomials  with 
positive  coefficients,  A(s)  and  D(s)  are  even 
functions  of  s and  of  degree  2n  and  2(n-l) 
respectively , 

B(s)  and  C(s)  are  odd  functions  of  s and  of  degree 
2n-l , 

(3)  Zeros  of  A(s) , [z  ] and  zeros  of  B(s) , [ z,  ] inter 


lace  on  the  jw  axis, 


bi 


zb  <|za 
D0  al 


< < • • •<  z. 


< Z. 


n-1 


n 


then 


(a)  zeros  of  A(s)  and  C(s)  interlace  on  the  jw  axis 

0 " zc  <lZa  Mzc  I < * * * < I 2C  I < I 2a  I and 
c0  al  C1  cn-l  an 


r 
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(b)  zeros  of  B(s)  and  D(s)  interlace  on  the  jw  axis 
0 - zw  <|zj  |<|zu  (<•••<  | Zj  | < | | and 


zb  <|zd  l<|zb  !<,,,<lzd 
b0  dl  °1  dn-l 


n-1 


(c)  zeros  of  C(s)  and  D(s)  interlace  on  the  ju>  axis 

°*z  <lzd  Nz,  l<,,,<lzd  M*c  I 

c0  dl  C1  dn-l  cn-l 

Proof  of  (a):  It  follows  from  (1)  that 

A(s)D(s)  - B (s)C  (s)  = 1 for  all  s, 
which  yields 

A(s)D(s)  - s2F(s)C(s)  = 1, 

where  sF(s)  = B(s)  and  sCT(s)  = C(s). 

„ , 2 
Substituting  z = s gives 

A(z)D(z)  - z F(z)C(z)  = 1 

where  zeros  of  A(z),  z and  F(z),  z . interlace  on  the 

ai  Di 

negative  Teal  axis  as  shown  in  Figure  2.1,  with 


(2-1) 


0>z  >z,  > 

ai  bi 


>z,  >z  , where  z * - z„ 
bn-l  an  ai  ai 


From  Equation  (2-1),  it  follows  that 

A(z  )D(z  ) - z 5(z  )C(z  ) = -z  F(z  )C(z  ) = 1 
i i i i i i i i 


and  similarly, 


-z  F(z  )C(z  ) 
ai+l  ai+l  ai  + l 


Therefore , 


&(*  )C(z  ) - za  F(z  )C(z 


ai*l  ai*l 


ai  + l 


(2-2) 


F(z  ) and  F(z  ) are  of  opposite  signs  because  of 
ai  ai*l 

the  zero  distributions  of  A(z)  and  F(z)  as  shown  in  Figure 

2.2.  Therefore,  from  Equation  (2-2),  C(z„  ) and  Cfz-  ) 

ai  ai+l 

are  also  of  opposite  signs.  This  implies  that  the  interval 


s 


(z  , z ) contains  an  odd  number  of  zeros  of  C(z).  There 
ai  ai  + l 

are  n-1  of  these  intervals  each  of  which  contains  one  and 

only  one  zero  of  C(z).  Therefore,  zeros  of  A(z)  and  C( z ) 

interlace  on  the  negative  real  axis  of  the  z-plane  with 

0 > z > z >..>  z > z„ 

al  C1  cn-l  an 

This  implies  that  zeros  of  A(s)  and  C(s)  interlace  on  the 
j bi  axis  with 

0<|za  !<|zc  l < * • * < I |<|z_  I 

al  C1  cn-l  an 

Since  s(T(s)  = C(s),  this  completes  the  proof  of  (a)  of 
Lemma  1.  Proof  of  (b)  and  (c)  follow  similarly. 

Q.E.D. 

Lemma  2,  3,  4,  and  5 are  summarized  in  Table  2.1 
without  proofs.  The  following  example  illustrates  the  zero 
interlacing  property  of  a matrix  which  satisfies  Lemma  1. 


Example  2.1.  Let 


Sn(s)  - 


6s4+6s2+l  6s3+3s 

C(s)  D(s) 


and  det  S fs)  * 1. 


As  can  be  seen  A(s)  and  B(s)  are  of  degree  four  and  three 

respectively.  The  zeros  of  A(s)  are  at  ±j  0.459,  ±j  0.888 

and  ther  zeros  of  B(s)  are  at  0,  ±j  0.707.  It  follows  from 

A(s)D(s)  - B(s)C(s)  3 1 that  C(s)  3 6s3  + 4s  and  D(s)  = 

2 

6s  ♦ 1 . The  zeros  of  C(s)  are  at  0,  ±j  0,816  and  the 
zeros  of  D(s)  are  at  ±j  0.409.  Indeed  the  zeros  of  A(s) 
and  B(s) , A(s)  and  C(s),  B(s)  and  D(s),  C(s)  and  D(s) 
alternate  on  the  j<*>  axis  as  stipulated  in  Lemma  1. 


4 


11 


12 


Sn(s) 


A(s)  B (s) 
C (s)  D(s) 


where 


A(s)  - a2ns2n*a2(n-l)sZln’1)*-"*  1 

B(s)  - blS 

C(s)  * c2n-ls2nl*c2n-3sZn^+‘'*+  cls 
D(s)  ‘ d2(n-l)s2(n'1)*d2(n-2)s2(n'2)*--'*  1 


Define  K (s)  * 
n 


1+s‘-lc  si 
n n n 


sc 


n 


Therefore  Kn’*(s)  » 


(2-4) 


— 

- 

- 

— 

— 

— 

1 

~s^n 

1 

0 

1 

-stn 

’scn 

l+s  Z£  c 
n n 

;scn 

1 

0 

1 

Let  Sn-l(s) 


A(s)  B(s) 
C(s)  D(s) 

A(s)  B (s) 
C(s)  D(s) 


* <*> 


-scn  1 


1 -S  L 


n 


(2-5) 


It  is  next  shown  that 


satisfies  conditions 


A(s)  B(s) 

|_C(s)  D(s)J 

(a) , (b) , (c) , (d)  and  each  of  its  entries  is  two  degrees 


lower  than  those  of  S fs)<‘s. 

n v 


i 


13 


Equation  (2-5)  can  be  written  as 


A (s)  B (s)  A(s)  - scnB(s)  B (s)  1 -s£n 


C (s ) D(s)  0 


C(s)  D(s)  C (s)  - scnD(s)  D (s ) 0 


A(s)  B(s)  1 -st 


(2-6) 


where  A(s)  = A(s)  - scnB(s) , 


(2-7) 


C(s)  = C(s)  - scnD(s) . 


(2-8) 


Equation  (2-7)  can  be  written  as 

■ <a2n  - Sb2n-l>s2n  * <a2n-2  ‘ c„b2„-3>s2n'2 


(a2n-4  • cnb2n-5)s2n'4  1 


(2-9) 


Choose  c. 


2n-l 


(2-10) 


This  causes  the  coefficient  of  in  Equation  (2-9)  to 
vanish.  Next  it  is  intended  to  show  that  the  coefficient  of 
s2n-2  Equation  (2-9)  is  positive,  i.e., 


a2n-2  ‘ cnb2n- 3 


a2n-2b2n-l  ' a2nb2n-3 


2n-l 


(2-11) 


It  follows  from  hypothesis  that 

A(s)  - a2ns2"  . a2n.2s2"'2  «•••*  a2s2  * 1 

’ b2n-ls2n'*  * b2n-3s2n'J  b3»5  * »ls 

’ s(b2a-ls2n'2  * b2n-3s2n"4 V2  * »1> 


s B(s) 


4 
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Substituting  z = sz  yields 


n 

n- 1 

2nz  a2n 

- 2Z  +’* 

n- 1 

. ,n-2 

T , z ♦ 

2n-l 

b,  ,z 
zn  - 3 

By  hypothesis,  it  follows  that  zeros  of  A(z)  , z&  and  B(z), 
zb  interlace  on  the  negative  real  axis,  as  shown  in 


Figure  2.1,  with 


0 > z„  > z,  > 
al  bl 


> z.  > z , where  zo  = -\z„ 
bn- 1 an  ai  ai 


From  the  zero  distributions  of  A(z)  and  B(z)  , it  follows 
that  the  sum  of  the  zeros  of  A(z)  is  less  than  the  sum  of 
the  zeros  of  B(z) . Therefore,  the  negative  sum  of  the 
zeros  of  A(z)  is  larger  than  the  negative  sum  of  the  zeros 
of  l(z) , 

a2n-  2 _ b2n-3 


2n 


> F 


2n-l 


(2-12) 


which  implies  ^2n-2°2n-\  ' a2nb2n-3  > °* 

?n  - A 

The  coefficient  of  s of  A(s)  is  positive  because 
the  negative  product  of  the  zeros  of  A(z)  is  larger  than 
the  negative  product  of  the  zeros  of  l3(z) . Similarly,  all 
other  terms  of  A(s)  are  positive.  Therefore  A(s)  is  an 
even  function  of  s and  of  degree  2(n-l). 

Next,  it  is  shown  that  C(s)  is  an  odd  function  of  s 
and  of  degree  2n-3.  Equation  (2-8)  can  be  written  as 

= (c2n-l'Cnd2n-2)s  + (-c2n-3‘cnd2n-4'1  s 


2n-5 


+ (c2n-5‘cnd2n-6)s  +. • •+  1. 


(2-13) 


4 


Since  A(s)D(s)  - B(s)C(s)  = 1,  the  coefficient  of 


s4n  2 must  be 


zero , i . e , , 
a2nd2n- 2 


b2n-lc2n-l 


0 


which  implies  that  cn 


2n 


'2n- 1 


2n- 1 
2n-  2 


(2-14) 


From  Equation  (2-14)  it  follows  that  the  coefficient  of 

s2nl  in  Equation  (2-13)  vanishes. 

The  coefficient  of  s2n  ^ of  C(s)  is 

j „ C2n- 1 ^ 

c2n-3  ' cn  2n-4  c2n-3  " d2n_2  * 2n-4 

c2n- 3d2n- 2 ' c2n-ld2n-4  „ n 

- — -j > u » 

d2n-2 

which  follows  from  the  zero  distributions  of  C(s)  and  D(s). 
Similarly,  all  other  terms  of  C(s)  are  positive.  Therefore 
C(s)  is  an  odd  function  of  s and  of  degree  2n-3. 

Next  it  is  shown  that 

B(s)  is  an  odd  function  of  s and  of  degree  2n-3. 

D(s)  is  an  even  function  of  s and  of  degree  2n-4. 

The  matrix 

of  A (s)  and  C(s)  interlace  on  the  ju>  axis  in  a pattern 
given  in  Lemma  2. 

It  follows  from  Equation  (2-6)  that 


A (s)  B(s) 
C(s)  D (s) 


now  satisfies  Lemma  2;  zeros 


Sn-1^ 


A(s)  B(s)  A (s)  B(s)  - sf  A(s) 


C (s)  D(s) J C (s)  D(s)  - s£nC(s) 


where  B(s)  = B(s)  - sfnA(s) 


D(s)  = D(s)  - s£nC(s) 


Equation  (2-15)  can  be  written  as 


(2-15) 


(2-16) 


B(s)  = (^ncnb2n.1-^na2n)s2n  1 ♦ (^ncnb2n-3+b2n-l'ina2n-2)sZn  1 


2n-  3 


+ (£ncnb2n-5+b2n-3'*na2n-4)s  +'"+  (bl'*n)s 


(2-17) 


The  coefficient  of  s2n  + 1 of  B(s)  is 


f_Cnh?n-1  ' f — a 9 n = Ln  K * , - L a~  =>  0. 

n n Zn  1 n Z n 2 n * 1 n 2n 


The  coefficient  of  " of  B(s)  is 


^ncnb2n- 3 + b2n-l  ‘ *na2n-2 


"n  |b2n^T  ' b2n'3  ' *2n 


-2  + b2n- 1 


# a2nb2n-3  ' a2n-2b2n-l  . L 
' » 2n-r 


Choosing  l = — c > 0, 

n a2n-2b2n-l  ' a2nb2n-3 


(2-18) 


(2-19) 


results  in  Equation  (2-18)  to  vanish.  Therefore,  the  degree 
of  B(s)  is  < 2n-3. 
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liquation  (2-16)  can  be  written  as 


D(s)  ‘ <Vnd2n-!-V2..1>s!"  * (Vnd2n-4*d2„-2- V2n-3>s 


2n-2 


♦ A + r.  C~  r)s^n'^  1 

n n 2n-6  2n-4  n 2n-5' 


(2-20) 


The  coefficient  of  s2n  of  D(s)  is 


Vnd2n-2  ' *nC2n-l  = ln  * d“7  * d2n-2  * *nc2n-l  = °‘ 


The  coefficient  of  s2n~2  of  D(s)  is 


£ c d,„  . + d, 
n n 2n - 4 2n  - 2 


^nc2n-3 


= l 


2n 


'2n- 1 


d2n-4  " c2n-3 


+ d 


2n-2 


2n- 1 


a2n-2b2n-l*  a2nb2n-3 


a2nd2n-4  ' b2n-lc2n-3 
b2n- 1 


+ d 


2n-  2 


b2n-l(a2nd2n-4 


b2n-lc2n-3)  + d2n-?(a 


2n- 2 ’•“2n- 2u2n-l  " a2nb2n-3) 


a2n-2b2n-l  ‘ a2nb2n-3 


(2-21) 


Now,  again  since  A(s)D(s)  - B(s)C(s)  = 1,  the  coeffi 
cient  of  s4n  4 must  vanish  which  implies  that 


a2nd2n-4  * b2n-lc2n-3  = b2n-3c2n-l  " a2n-2d2n-2•  (2-22) 


Substituting  Equations  (2-14)  and  (2-22)  into  (2-21)  causes 
the  coefficient  of  s2n'2  of  D(s)  to  vanish.  Therefore,  the 
degree  of  D(s)  is  < 2n-4. 
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Now,  consider  the  matrix  Sn  ^(s)  ■ 


A(s)  B(s) 
C(s)  D(s) 


where 


A(s)  is  an  even  function  of  s and  of  degree  2n-2 
B(s)  is  an  odd  function  of  s and  of  degree  £ 


Ms) 

is 

an 

B(s) 

is 

an 

2n- 

■3 

C(s) 

is 

an 

D(s) 

is 

an 

2n- 

-4. 

C(s)  is  an  odd  function  of  s and  of  degree  2n-3 
D(s)  is  an  even  function  of  s and  of  degree 

Since 

A(s)D(s)  - B(s)C  (s)  * 1, 

A (s)D(s)  - s28(s)C(s)  - 1, 

where  both  B(s)  and  C(s)  are  even  functions  of  s. 

2 

Let  z = s . Then, 

A(z)D(z)  - z B(z)C(z)  * 1, 

and 


zfl  B (z  )C (z  ) = z B(z  )C(z  ).  (2-23) 

ai  ai  ai  ai+l  ai+l  ai+l 


C(z  ) and  C(z  ) are  of  opposite  signs  because  of 
ai  ai+l 

the  zero  distributions  of  A(z)  and  C(z).  Therefore, 

B(z  ) and  B(z  ) are  also  of  opposite  signs,  This 
ai  ai+l 

implies  that  the  interval  ( z , z ) contains  an  odd 

ai  ai+l 

number  of  zeros  of  B(z).  There  are  n-2  of  these  intervals 
each  of  which  contains  one  and  only  one  zero  of  B(z). 
Therefore  B(z)  is  of  degree  >_  n-2  and  B(s)  is  of  degree 
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\ 

i 

I 


i 


i 


> 2n-4  and  B(s)  * slf(s)  is  of  degree  >_  2n-3.  Since  the 
degree  of  B(s)  has  been  previously  shown  to  be  <_  2n-3, 
therefore  B(s)  is  of  degree  2n-3.  Similarly,  D(s)  is  of 
degree  2n-4.  This  proves  the  Sn_^(s).  This  decomposition 

process  is  continued  until  n=l.  This  results  in 


Define 


Sj  (s) 


v2*1 


C1S 


K.  (s) 


l+s^1c1 


SC, 


where 


Cj  and  Tj 


b 


1* 


The  product 
SjCsJK^Cs) 


a2s  +1  b^s 


cls 


(a2-b1c1)s^l 


(c1-c1)s 


1 -ll  s 

-C^S  ^^CjS^+1 

(bl£lcl'a2'ei^s3+(bl'£l)s 
(f1c1-c^1)s2  + l 


Therefore,  the  above  process  results  in 

V»>  - K1(sV1C2(s)...1Cii(s) 

as  described  in  the  theorem.  This  completes  the  proof  of 
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the  sufficient  part  of  the  theorem. 

Necessity:  It  is  intended  to  prove  by  induction  that  if 

there  is  Kj (s) ‘K2 (s) • • *K  (s)  where  each  KjCs)  is  described 


by  Equation  (2-3),  then 


Sw(») 


A(s)  B(s) 
C (s)  D(s) 


possesses  the  properties  (a) , (b) , (c) , and  (d)  of  the 
theorem. 

Obviously  Sj^s),  S2(s)  •••  satisfy  conditions  (a),  (b)  , 

(c)  , and  (d) . If  Sn(s)  has  properties  (a),  (b)  , (c)  , and 

(d)  , then 


Sn(s>Ws> 


A(s)  B (s)  1+s‘l^.c 


C(s)  D(s) 


'n+1  n+1  n+1 


A (s)  BCs)  1 0 

C (s)  DCs)  scn+1  1 


_ACs)  BCs)  1 T 1 s£n  + 1 1 0 

C (s)  D(s)  J 0 1 J scn+1  1 


(2-24) 


Then,  B(s)  * s£n+1A(s)  ♦ B(s)  and  of  degree  2n+l, 
D(s)  * s£n+1C(s)  + D(s)  and  of  degree  2n, 


Since  det 


A(s)  B(s) 
_C(s)  D(s) 


* 1 and  zeros  of  A(s)  and  C(s)  inter- 


lace on  the  jw  axis,  it  follows  that,  by  Lemma  4,  zeros  of 
B(s)  and  D(s)  interlace  on  the  ju  axis. 


d 
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Equation  (2-24)  is  now  written  as 


S„(s)Kn+1(s) 


n 


A(s)  B (s) 
C(s)  D(s) 


f 


where  A(s)  * A(s)  + scn+1B(s)  and  is  of  degree  2(n+l), 
C(s)  ■=  C(s)  + scn+jD(s)  and  is  of  degree  2n+l. 


Since  det 


A (s ) 
C(s) 


B(s) 

D(s) 


1 and  zeros  of  B(s)  and  D(s)  inter- 


lace on  the  ju>  axis,  it  follows  that,  by  Lemma  5,  zeros  of 
A(s)  and  B(s)  , A(s)  and  C(s),  C(s)  and  D(s)  interlace  on 
the  jw  axis.  This  shows  that  the  matrix  Sn+1(s)  also  has 

properties  (a),  (b) , (c) , and  (d)  of  the  theorem. 


Q.E.D. 


Corollaries  to  this  theorem  are  summarized  in  Table 
2.2  without  proofs.  Theorem  2.1  and  its  corollaries  for- 
mulate the  method  of  decomposition  from  the  right.  The 
following  theorem  shows  decomposition  from  the  left. 


Theorem  2.2  Let  the  matrix  S (s)  * 

n 


A (s ) B (s) 
C (s)  D (s) 


Sn(s)  - Kn(s)Kn.1Cs),,*K2K1 


where 


K.(s) 


1+s^.^  sti 


sci 


, l^,  ci  > 0 for  i - l,2,***n 


Table  2.2  Sununai^  of  Decomposition  Corollaries 


Zeros  of  C(s)  and  D(s)  interlace  Zeros  of  A(s)  and  B(s)  interlace 
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if  and  only  if 

(a)  Det  Sn(s)  = 1, 

(b)  A(s),  B(s),  C(s)  and  D(s)  are  polynomials  with  positive 
coefficients  A(s)  and  D(s)  are  even  functions  of  s 

and  of  degree  2n  and  2(n-l)  respectively, 

B(s)  and  C(s)  are  odd  functions  of  s and  of  degree 
2n- 1 , 

(c)  Zeros  of  A(s),  [z  ] and  zeros  of  C(s),  [z  ] interlace 

ai  ci 

on  the  joi  axis  with 


0 


lz. 


'n-1 


(d)  A (0)  = 1. 

Proof  of  this  theorem  is  similar  to  that  of  Theorem  2.1  and 
therefore  omitted. 


Q.E  .1), 

Consider  a passive,  reciprocal  2-port  network  with  the 
polarities  of  Ij,  I2,  , and  as  shown  in  Figure  2.3. 

The  chain  matrix  expression  of  this  2-port  is  then  given  as 


A B 
C D 


I. 


If  this  2-port  is  of  simple  ladder  structure  of  Figure  2.4 
then  the  chain  matrix  of  this  ladder  is 


1 z 

1 0 

1 +YZ  Z 

0 1 

Y 1 

Y 1 

■ 


If  all  series  arms  and  shunt  arms  are  inductive  and 
capacitive  elements,  respectively,  then  Z(s)  = Is, 

Y(s)  - cs,  and  therefore  matrix  in  the  foregoing 

theorem  is  realizable  by  this  ladder.  It  follows  that  a 
cascade  of  n of  these  ladders  as  shown  in  Figure  2.5  has 
an  overall  chain  matrix  of  K^«K2*  »»*K  ■ S which  is  the 

overall  matrix  stated  in  the  theorem. 

Section  2,4  The  Decomposition  Algorithm 

The  algorithm  for  determining  the  values  of  ladder 
elements  will  now  be  established.  The  decomposition  theorem 
is  used  to  develop  a computer -supported  procedure  based  on 
the  two  known  entries  in  the  overall  chain  matrix. 

Ms)  - a2ns2n  + a2(n.1)sZ(n-1)  +•••♦  a2s2  ♦ 1 

B(s)  “ b2n-ls2n  1 + b2n-3s2n'3  +**,+  b3s3  + bls 
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1 


A Routh 

array  is 

formed ; 

a2n 

b2n-l 

a2n 

®2Cn-l) ’ * * *2  1 

s2n 

b2n-l 

b2n-3  * bl 

2n-l 

s 

a2 (n-1) 

°2n-l 

°2(n-l) 

Y2n-3 

a2  (n-1) 

a2 (n-2)  " ’ a0 

s2(n-l) 

Y2n-3 

• 

T2n-s’  bl 

2n-3 

s 

• 

• 

t 

• 

y 

• 

where  each  element  starting  from  the  second  column  to  the 
right  is  evaluated  by  Routh  criterion  method,  i.e.f 


b2n-la2(n-l)  ' a2nb2n-3 
°2  (n-1)  = ~ ^Tri-T” 


b2n-la2 (n-2)  ~ a2nb2n-5 
a2 (n- 2)  “ b2n-l 

4 a2(n-l)b2n-3  ' b2n-la2(n-2) 
Yzn*3  * a2 (n- 1) 


, a2(n-l)b2n-5  ' b2n-la2(n-3)  etc 
Y2n'5  °2 (n-1) 


1 


■>» 
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In  the  first  column,  (k,l)  entry  gives  the  value  of 
ly  if  k is  odd  and  the  value  of  if  k is  even.  All  the 

elements  in  the  array  and  the  values  of  I's  and  c’s  can  be 
computed  by  a simple  digital  computer  program. 

The  above  algorithm  applies  to  the  circuit  as  shown  in 
Figure  2.6  where  the  number  of  reactive  elements  is  even 
(Theorem  2.1).  If  the  number  of  reactive  elements  is  odd  as 
shown  in  Figure  2.7  (Corollary  2,2),  then  the  algorithm 
should  be  modified  as  follows: 

(1)  In  the  Routh  array,  the  coefficients  of  B(s) 
should  be  located  in  the  first  row,  and  those  of  A(s)  in 
the  second  row. 

(2)  In  the  first  column,  (k , 1 ) entry  gives  the  values 
of  if  k is  odd  and  the  values  of  if  k is  even. 

In  both  Figure  2.6  and  Figure  2.7,  m denotes  total 
number  of  elements. 

Example  2.2 

A(s)  - (s2  + l)  (10s2*l) (20s2  + 1) (36s2  + 1) 

- 7200s8+8460s6+1346s4+67s2+1 

B(s)  - (6s2+1)(12s2+1)(25s2+1)s 

- 1800s7+522s5+43s3+s 
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Cj  = 4 

7200 

8460 

1346 

67 

1 

s8 

1 

t2  - 0.28 

1800 

522 

43 

1 

0 

S7 

C3  - 33.47 

6372 

1174 

63 

1 

0 

S6 

l 4 3 0.58 

190.36 

25.20 

0.72 

0 

0 

s5 

c5  = 120.53 

330.36 

38,98 

1 

0 

0 

s4 

*6  " °*12 

2,74 

0.14 

0 

0 

0 

s3 

c?  = 1335.54 

21.95 

1 

0 

0 

0 

s2 

l8  = 0.016 

0.016 

0 

0 

0 

0 

s1 

1 

0 

0 

0 

0 

S° 

The  completed 

circuit 

is  shown 

in  Figure  2. 

8. 

Example  2.3 

(Corollary  2,2) 

Given: 

A(s)  - 240s6  + 234s4  + 

48s2  + 

1 

B(s)  =■  1440s7  + 1452s5 

+ 330s3 

+ 

12s 

l1  = 6 1440  1452 

330 

12 

s7 

c2  3 5 240 

234 

48 

1 

s6 

(s* 

II 

ro 

00 

42 

6 

s5 

c4  3 4 24 

18 

1 

s4 

*5  " 3 6 

4 

s3 

' 

c6  = 2 2 

1 

s2 

l7  3 1 1 S1 


* 

The  completed  circuit  is  shown  in  Figure  2,9, 


CHAPTER  3 


I 

1 


APPLICATION  OF  CHAIN  MATRIX  DECOMPOSITION  TECHNIQUES 

Section  3.1  Introduction 

In  this  chapter  the  techniques  which  were  developed 
earlier  will  now  be  used  in  several  practical  situations 
which  are  of  interest.  Fortunately,  the  literature  [25]  is 
rich  in  characteristics  for  filter  networks  in  wide  variety 
and  degrees  of  sophistication.  For  example  such  filters  as 
Butterworth,  Chebyshev,  and  Bessel  are  well  known  and  com- 
pletely characterized.  These  will  now  be  synthesized 
through  use  of  the  decomposition  of  chain  matrices.  The 
elements  of  these  matrices  are  predetermined  by  the  ortho- 
gonal polynomials  used  to  approximate  the  ideal  filter 
characteristics.  The  chain-matrix  method  is  shown  to  be 
simple  and  straightforward.  Results  of  the  above  filters 
are  compiled  into  a table  for  both  single  and  double-ter- 
minated cases. 
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Section  3.2  Synthesis  of  Single-Terminated  Ladder  Network 
by  Chain  Matrix  Decomposition 
An  ideal  low-pass  filter  characteristic  of  Figure  3.1 
can  be  approximated  in  various  ways,  such  as 


V2(s) 


V^s) 


1 + U) 


7n 


s = j OJ 


(3-1) 


f 


A 


34 


Vj(s) 


(3-2) 


2 

Vx(s) 

S = ju) 


(3-3) 


where  Vj(s)  and  V2(s)  are  the  input  and  output  voltages  of 

the  filter  networks  respectively. 

As  is  well  known.  Equations  (3-1),  (3-2),  and  (3-3) 
result  in  Butterworth,  Chebyshev,  and  Bessel  polynomials 
and  hence  the  corresponding  filters.  Synthesis  of  these 
filters  have  been  extensively  investigated  and  widely 
documented.  They  are  discussed  and  illustrated  in  most  of 
the  textbooks  on  network  synthesis.  In  this  section, 
synthesis  of  these  filter  networks  are  accomplished  through 
the  application  of  the  decomposition  of  chain  parameter 
matrices  whose  elements  are  predetermined  by  the  orthogonal 
polynomials  used  to  approximate  the  ideal  characteristics. 
These  polynomials  are  well  known  in  the  literature  and  can  be 
found  in  many  mathematical  handbooks. 

Consider  the  circuit  shown  in  Figure  3.2  in  which  m 
number  of  reactive  elements  are  connected  to  a voltage  source 
and  terminated  in  a 1-Jl  resistor.  The  chain  parameter  matrix 
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Since  V2(s)  = I2(s),  it  follows  from  Equation  (3-4)  that 

V2(s)  j 


V^s)  A(s)+B(s) 


(3-5) 


Therefore , 


V2(s) 

2 

s 

1 

1 

V j ( s ) 

s*  ju 

A(s)+B(s) 

A(-s)+B(-s) 

s»ju 


= 1 1 

V5)  B«<-» 

where  Bm(s)  are  orthogonal  polynomials  of  mth  order  used 

for  the  approximation  of  the  ideal  filter  characteristics. 
It  follows  from  Equation  (3-6)  that 

A - -L_. 

A (s) +B(s)  Bm(s) 

Since  A(s)  and  B(s)  are  the  elements  of  the  chain  matrix  of 

Theorem  2.1,  A(s)  is  an  even  polynomial  and  B(s)  is  an  odd 

polynomial  and  hence  they  are  uniquely  determined.  A(s)  * 

Ev [ B ( s ) ] , B(s)  = Od[B_(s)].  In  other  words,  if 
m hi 

Bm(s)  = sm+am.1Sm'1+am.2sm'2+...+a1s+a0  (3-7) 

i i 

where  m is  even,  then 


(3-6) 

S = ju> 


I 

t 


A(s) 


sm+a  7sm-Z+. 
m-  L 


+a2s  +a0 
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(3-8) 


sm‘3*. 


+ajS  +ajS 


(3-9) 


The  following  examples  demonstrate  the  application  of  de- 
composition method  to  filter  designs. 

Example  3.1.  (Application  of  Theorem  2.1) 

Given  6th  order  Bessel  polynomial 

B6(S)  = s6+21s5+210s4+1,260s3+4,725s2+10,395s+10,395 


Solution : 

A (s)  = Ev[B6(s)]  = s6+210s4+4,725s2+10,395 
B (s)  * 0d[B6(s)]  = 21s5+1,260s3+10,395s. 

A Routh  array  is  formed  and  circuit  elements  computed  in 
accordance  with  the  deomposition  algorithm  developed  in 
Chapter  2: 


c6=0. 048 

1 
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4725 

10395 

S° 

c n 

II 

O 

21 

1260 

10395 

0 

s5 

c^=0 . 225 
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4230 

10395 

0 

s4 

f 2=0 . 301 

667.8 

8939.7 

0 

0 

s3 

c2=0 . 382 

2221.98 

10395 

0 

0 

2 

s 

f 1=0. 559 

5815.56 

0 

0 

0 

s1 

10395 

0 

0 

0 

s° 

i 


; 
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The  completed  6th  order  Bessel  filter  circuit  is  shown  in 
Figure  3.3. 

Example  3.2.  (Application  of  Corollary  2.2) 

Given  5th  order  Butterworth  polynomial 


B5(s)  - 

s5+3.236s4+5 

•236s3+5.236s2+3.236s+l 

Solution : 

A(s) 

- nv[B5(s)] 

= 3.236s4+5.236s2+1 

B(s) 

= 0d[B5(s)] 

- s5+5.236s3+3.236s. 

£s=0.309 

1 

5.236  3.236  s5 

c4=0.894 

3.236 

5.236  1 s4 

£ j=1 . 382 

3.618 

2.927  s3 

c2 =1.694 

2.618 

1 s2 

£1=1.545 

1.545 

s1 

r 

1 

s° 

[ 

The  completed  5th  order  Butterworth  filter  circuit  is 

1 . 

shown  in  Figure 

Tables  3.1 

3.4. 

, 3.2,  and  3 

.3  show  the  normalized  element 

values  for  a single- terminated  Butterworth,  Chebyshev,  and 
Bessel  filters  respectively. 


Table  3.1  Normalized  Element  Values  for 


Table  3.3  Normalized  Element  Values  for 
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Section  3.3  Synthesis  of  Double-Terminated  Ladder  Network 
by  Chain  Matrix  Decomposition 
Consider  the  circuit  shown  in  Figure  3.5.  This  lossless 
ladder  network  is  described  by  chain  parameters  A(s),  B(s), 
C(s),  and  D(s)  and  it  is  double- terminated  by  one-ohm  re- 
sistors. Its  chain  matrix  is 


V(s) 

I^s) 


A(s) 

C(s) 


V s) 

vs)_ 

(3-10) 


1 

t/> 

w 

rH 

‘1 

A (s)  B(s) 

V2(s) 

i 

H 

o 

1 

C (s)  D(s) 

_I2(S). 

It  follows  from  Equation  (3-11)  that 

Vj (s)  = [A(s)+C(s)]V2(s)+[B(s)+D(s)]I2(s). 

Since  V2(s)  = I^fs),  Equation  (3-12)  results  in 

v2(s)  l 

vprr  " A(s)TB(S)+cTi)+DrsT- 


(3-11) 


(3-12) 


(3-13) 


The  voltage  transfer  function  of  the  ladder  network  of 
Figure  3.5  can  be  written  as 

V2  (s ) H 

VJTsT  ’ Bj^TsT’  (3'14 

whore  B (s)  is  the  normalized  orthogonal  polynomial  of  mth 
m 

order  that  is  used  to  approximate  ideal  low-pass  filter 
characteristics  and  H is  the  constant  multiplier.  At  s = 0, 
the  network  is  purely  resistive  and  therefore. 
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V2(s) 


Vj(s) 


s = 0 


H 

Bmt°) 


1 

1 


Since  normalized  BfO)  ■ 1,  the  constant  multiplier  H is 

m 

one  half.  Thus 


Vs)  . 1 

Vj(s)  A(s)+B(s)+C(s)+D(s) 


2B  (s) 
m 


(3-15) 


In  Chapter  2,  it  is  shown  that  for  an  LC  ladder  network, 
A(s)  and  D(s)  are  even  functions  of  s and  B(s)  and  C(s)  are 
odd  functions  of  s.  Therefore,  it  follows  from  Equation 
(3-15)  that 

A(s)+D(s)  - Ev(2Bm(s)]  (3-16) 

B(s)+C(s)  - 0d[2Bm(s)] . (3-17) 


The  objective  of  this  section  is  to  apply  the  chain 
matrix  decomposition  method  to  synthesize  double- terminated 
Butterworth  filters.  The  terminations  are  assumed  to  be 
one  ohm  as  shown  in  Figure  3.5.  It  has  been  shown  (25]  that 
these  filters  are  symmetrical  if  m is  odd.  Furthermore,  if 
m is  even,  they  are  antimetrical . Therefore, 

A(s)  = D(s)  for  m odd,  and  (3-18) 

B(s)  ■ C(s)  for  m even.  (3-19) 

Equation  (3-15)  to  (3-19)  are  useful  for  synthesizing 
double-terminated  filters  by  chain  matrix  decomposition. 
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The  following  example  illustrates  the  procedure. 

Example  3.3:  Synthesize  the  double-terminated  6th  order 
Butterworth  filter  given  the  Butterworth  polynomial. 

B6(s)  * s6+3.8637s5+7.4641s4+9.1416s3+ 
7.4641s2+3.8637s+1 

Solution:  From  Equations  (3-16)  and  (3-17),  it  follows 
that 

A(s)+D(s)  = Ev[2B6(s)]  = 2(s6+7.4641s4+7.4641sZ+1)  (3-20) 

and 

B(s)+C(s)  = Od[2B6(s)]  = 2(3. 8637s5+9.1416s3+3. 8637s) . 


Since  m=6 , it  follows  from  Equation  (3-19)  that 

B (s)  = C(s)  * 3.8637s5+9.1416s3+3.8637s.  (3-21) 

For  passive  networks,  A(s)D(s) -B(s)C(s)  = 1,  thus, 

A(s)D(s)  - B (s)C (s) +1 

- (3.8637s5+9.1416s3+3.8637s)2+l. 

But , 

(A(s)-D(s)]2  = [A(s)+D(s)]2-4A(s)D(s) 

* 4(s6+7.4641s4+7.4641s2+1)2- 
4(3.8637s5+9.1416s3+3.8637s)2-4 


Therefore , 
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A ( s ) -D(s)  =»  2s6.  (3-22) 

Adding  Equations  (3-20)  and  (3-22)  results  in 

A(s)  « 2s6+7.4641s4+7.4641s2+l.  (3-23) 

The  chain  parameters  A(s)  [Equation  (3-23)]  and  B(s) 
[Equation  (3-21)]  can  now  be  used  to  realize  the  filter. 

The  computation  of  element  values  is  shown  in  Table  3.4. 

The  corresponding  filter  network  is  shown  in  Figure  3.6. 
Table  3.5  shows  the  normalized  element  values  for  a double- 
terminated  Butterworth  filter. 
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Table  3.4  Computation  of  Element  Values 
of  Example  3.5 


0.5176 

1.4142 

1.9319 

1.9319 

1.4142 

0.5176 


2 

7.4641 

7.4641  1 

3.8637 

9.1416 

3.8637 

2.7321 

5.4641 

1 

1.4142 

2.4495 

0.7320 

1 

0.5176 

Figure  3.4  5th  Order  Butterworth  Filter 


A,B,C,D 


Figure  3.5  General  Form  of  Double-Terminated  Filter 


Figure  3.6  6th  Order  Butterworth  Filter 


CHAPTER  4 


OPTIMAL  SYNTHESIS  OF  INHOMOGENEOUS  LADDER 
NETWORK  USING  CHAIN  MATRIX 

Section  4.1  Introduction 

The  synthesis  of  inhomogeneous  ladder  networks  using 
chain  matrix  techniques  has  not  been  investigated  to  any 
significant  extent.  Very  little  material  has  appeared  in 
the  literature  on  this  subject.  A recent  paper  by  Lee  and 
Brown  sheds  some  light  on  the  pole-zero  distribution  pattern 
of  the  immittance  function  for  inhomogeneous  ladder  networks. 
Much  more  work  is  necessary  in  order  to  clearly  characterize 
such  systems.  In  this  chapter  several  synthesis  procedures 
of  optimal  inhomogeneous  ladder  networks  are  developed 
through  use  of  the  pole- zero  distribution  pattern  developed 
by  Lee  and  Brown  and  a simple  complex  variable  transforma- 
tion. 

Section  4.2  Optimal  Synthesis  of  Inhomogeneous  Ladder  Networks 
Fundamental  definitions  concerning  inhomogeneous  ladder 
networks  and  optimal  synthesis  are  given  below  in  order  to 
facilitate  proofs  of  Lemmas  which  will  establish  the 
synthesis  procedures. 

Definition  4.1:  A ladder  network  is  said  to  be  inhomogenous 
if  all  its  series  arm  impedances  are  f^z(s)  and  all  its 

shunt  arm  admittances  are  g^y(s),  where  i * 1,  2,  •••,  n and 
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z(s)  f y(s).  Such  a network  appears  in  Figure  4.1- 

Definition  4.2  An  arc  is  denoted  by  fG(s),  k] , if  it 
satisfies  the  root  locus  equation, 


1 + — — = 0, 
G(s) 


for  k > 0 in  the  s-plane. 


Definition  4.3  Suppose, 


P(s)  - k TT  [G (s)  +p . ] 

" i=i  1 


N(or  N-L) 

Q(s)  = k TT  [G (s) +q. ] 
q i-1  1 


If  0 < pj  < q^  < P^+1  for  all  i = 1,  2,  N-l,  then  zeros 

of  F(s)  and  Q(s)  are  said  to  interlace  with  respect  to 
[G(s) , k] . [Observation  is  made  here  that  the  zeros  of 
P(p)  and  Q(p)  interlace  on  the  negative  p = G(s)  axis  of  the 
plane. ] 

Definition  4.4  Synthesis  of  an  inhomogeneous  ladder  network 
from  a given  transfer  function  is  said  to  be  optimal  if 
(rf^fEgj)  is  minimized.  The  quantities  f^  g^^  are  the 

impedance  and  admittance  level  for  the  ith  section  of  the 
network  of  Figure  4.1. 


A(s)  B(s)  N 

Let  S (s)  * * TT  K.  be  the  chain  matrix  of 

n C (s)  D(s)  i-1  1 


Figure  4.1. 


If 


K. 


1+figiy(s)z(s)  £,.(.) 

giy(s)  i 


n (s ) n.(s) 

where  z(s)  = d ^y,  y(s)  * d ^ sy , n *=  degree  of  [na(s)nb(s)] 
a b 


or  degree  of  [da(s)db(s)]  whichever  is  greater.  Then 


Nn 

TT  (s+zfl  ) 
A(s)  - Ka  — ^ 


[da(s)db(s)]‘ 


(4-1) 


A(s) 


[s |na(s)nb(s)  « 0] 


z_1(s)B(s)  = Kb  lBl 


(N-l)n 

IT  (s+zK  ) 


[da(s)db(s)] 


nn:' 


(4-2) 


z* 1 (s) B(s) 


[s|na(s)nb(s)  « 0] 


*fi  * fT 


■ 
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y'A(s)C(s)  = Kc 


(N-l)n 
IT  (s  + z ) 
i=»l  ci 

[da(s)db(s)]N’1 


(4-3) 


y_1(s)C(s) 


[s |na(s)nb(s)  = 0] 


D(s) 


(N-l)n 

ff  (s+zd  ) 

Cd  [da(s)db(s)]ff"T* 


(4-4) 


D(s) 


(s |na(s)nb(s)  - 0] 


where  zeros  -zfl  and  -zb  , -zfl  and  -zQ  , -zb  and  -z^  , -zc 

and  -z,  , interlace  with  respect  to  [y(s)z(s),  k] . There- 
ai 

fore,  it  follows  from  observation  that  there  exists  a trans- 
formation of  y(s)z(s)  = p such  that 


A (p)  - TT  (1+£r>’ 

i-1  1 


N-l 


(z^BKp)  - fT  TT  (!♦£-), 

i-1  i 


1 


4 


[y'lc](p)  • eT  iY1(1’n) 


N-l  n 

D(p)  3 TT  (l+r-) » 
i-1  i 


where 


0 < < 6^  < ai+l» 

0 < oi  < < ai+l» 

0 < ei  < 6i  < ®i+l» 

0 < Yi  < fii  K Yi+1‘ 

If  an  inhomogeneous  ladder  network  ends  in  a series 
impedance  as  shown  in  Figure  4.2,  then  its  chain  parameters 
have  the  following  forms: 


A(S)  - Ka  — 


Nn 

TT  (s+z  ) 

i*l  ai 


(4-5) 


[da(s)db(s)1 


A(s) 


■ 1. 


ts|n  (s)nb(s)  - 0] 


z *^B (s)  - K 


Nn 

*TT  (s+zb  ) 
i-1  i 

b [da(s)db(s)]1 


(4-6) 


z* 1 (s)B(s) 


tft  - £t* 


[s|na(s)nb(s)  = 0) 
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y'1  (s)c(s) 


(N-l)n 

IT  Cs+z  ) 

i = l Ci 

[da(s)db(s)]N-: 


(4-7) 


y_1(s)C(s) 


fs|na(s)nb(s)  = 0] 


8 j 8j  i 


1)(S)  - K„ 


Nn 

.TT  (s+z.  ) 


fda(s^db(s^ 


(4-8) 


D(s)  = 1, 

[s|na(s)nb(s)  = 0] 


where  zeros  -z  and  -zb  , -z  and  -z  , -z,  and  -z,  , 
i Di  ai  ci  bi  di 

-zc  and  -Zj  , interlace  with  respect  to  [y (s) z (s)  ,k] . 

Therefore,  there  exists  a transformation  of  y(s)z(s)  = p 
such  that 

N 

A(P)  = TT  (l+-f), 
i = l °i 


u bkp)  = f r r d+#-), 

1 i=l  3i 


[y  lC]  (p)  = g TT  (l+£-)f 
1 i = l ' i 


D(P)  - TT  (l+£.)f 
i = l 6i 


where 


0 < «i  < 6i  < ai+l * 

0 < “i  < yt  < “i*i> 

0 ‘ «i  < 5i  ‘ <W 

0 < Ifj  < {i  < ri*l- 


Lemma  4.1:  Let  A(s),  B(s),  C(s),  and  D(s)  be  the  chain  para- 
meters of  any  2-port  network.  If 


A(s)  « K, 


Nn 

TT  (s*za  ) 
i = l *i 

[g(s)]N 


g1(s) 

fx(s) 


A(s)  = 1, 

ls|f(s)  = 01 


Nn 

TT  (s  + z.  ) 
i = l Di 


b 18(3)]“  • 


gx(s) 

^(s) 


f-r  » 


[s|f(s)  = 0] 


wheTe  fj(s)-f2(s)  * f(s),  g1(s)-g2(s)  = g(s),  -z& 
interlace  with  respect  to  1 M » tlien 


and  -z. 
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fx(s) 

gi  (s) 


and 


N 


C(s)  = pA  (p) 


i = l aiA^(-ai)[z  B]  (-cu)  (p+o^) 


D_f(s) 
P g(s) 


(4-9) 


-1, 


D(s)  = LL_AIi£l 


1 - f t 21 


tt 


i=l  Bi [ z iB]'(-6i)A(-Bi)(p+Bi) 


p-f(s) 

p jU) 


l - 10) 


where  -a^  and  -B^  are  the  zeros  of  A(p)  and  [z  *B] (p)  in 
p-plane. 

Proof  of  Lemma  4.1:  It  follows  from  the  hypothesis,  there 
exists  a transformation  p = y(s)z(s),  such  that 

N 

A (p)  - TT 

i = l ai 


N 


[ z*  JB] (p)  = fT  (1  + BL)’ 


i-1  wi 


where  0 < < B^  < ai  + j. 

The  open  circuit  input  admittance  is  given  by 


1 . [zC](p) 

pzu  p 

■ *1 

v 


4 


where 


57  hi 


i-1  p+V 


(4-11) 


, , [ZC](D) 

Ai  P'TpT 


px-Oi 


(2C](-a.) 


tzCj  (p)  x J ^c](~ai) 
A(p)  Pixj  -aiA'*  C - « ± ) (p  + ai 


Since  A(p)D(p)- [z’ AB] (p) [zC] (p)  - 1, 


then 


[zC] (-a.) 


[z’AB](-a.) 


Thus  Equation  (4-11)  becomes 


[zC] (p)  _ _ ? 1 

A'fn  l — P 1 n 

i=l  a.A^-a-Hz  AB]  (-<*.)(?♦<*.) 


N 1 

[zCJ(p)  = pA(p)  L 

i*l  a.A'(-a.)[z  1B](-a.)(p+ai) 


This  completes  the  proof  of  Equation  (4-9). 

To  prove  Equation  (4-10),  consider  the  short-circuit 


4 


input  admittance 


Y±1  . D (£) 

P pU'^BKp) 


. N l)(-B.j 

4-  y, 1 

fTP  i = l (p+0i) 


rTP  i = l + 


Rearranging  the  expression  for  yields: 


_ P(P) = 1 . p J i 

(7^B](p)  ’T  i = l BiA(-Bi)[z*TB]"(-B.)(p+8.) 


It  follows  therefore  that 


Dtp)  - JilVpiL,  ; — T e 

T i = l B.[z'iB]n-Bi)A(-Bi)(p+Bi) 


This  completes  the  proof  of  Equation  (4-10). 


Q.E.D. 


Lemma  4.2:  Let  A(s),  B(s),  C(s),  and  D(s)  be  the  chain 
parameters  of  the  network  of  Figure  4.2.  If 


4 
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Proof  of  Lemma  4.2:  Again  by  observation,  as  in  the  proof 
of  Lemma  4.1,  there  exists  a transformation  p - y(s)z(s), 
and  A(s ) and  z_1(s)B(s)  become 


A(p)  =»  TT  (l*-f-), 

i=l  °i 


[z_1B](p)  = fT  TT  (1+|-), 


T 11  ^ x a- 

1 i = l ei 


where  0 < a.  < 6^^  < aj  + j 


The  open-circuit  input  admittance  is  given  by 


1 

pznCpJ 


[zCj  (P) 

pA(p) 

N-l 

*Tpi71(1*n') 

R n 

P 11 

i = l “i 


. , p+a . ’ 

1 = 1 v l 


(4-14) 


where 


(zC](p) 

pA'TF)" 


p*  - a 


[zC] (*ai) 


-aiA-(-ai, 


P 


1 
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•••  ™ ’ W 


N | zC ) ( - a j ) 

l’i=i  r^Xp^sp' 


Since  A(p)D(p)  - [ z" AB] (p) [zC] (p)  = 1, 


then 


[zC](-a.) 


[z'1B](-ai) 


" 

A(p)  i = l c^A'f-a.) 


[z'^B]  (-a.)  Cp  + t*i) 


From  Equation  (4-14)  it  follows  that 


iis  W 


i«l  a^A'C-a.) [z'1B](-oi) 


(4-1S) 


This  completes  the  proof  of  Equation  (4-12). 

To  prove  Equation  (4-13),  consider  the  following 
fraction  expansion: 


pA(p)[z*1B](p) 


p * (l+-E-)fT  v (1+tP-) 
i=l  ai  ‘i=i  Bi 


, N K-  N K'a 

^ * ili  Pr»T  4 ill  l 
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11m  1 j 1 1 

p >0  " A(n)[z'TBJ  (p)  ^ 


N 

Z 


1 


i = l a.Va.Hz'^K'a.) 


N 

>: 


i=l  3iZA(-6i)fz"1Bl'(-Ri) 


N 


gj  - 2 j TT 

i=l  Bi/A(-6i) [z  IB]'(-Bi) 


On  the  other  hand: 


(4-16) 


1 


1 


A(p)|z‘JB|  (p)  i,(l+J?  )rTi,(]  + IL-) 


(l+ajp* • • • ) (fT+b1p+ • • • ) 


fT+^alfT+bl^p+ ^albl+’ * ' )p  +' 


1 1 °i 

t—  - T-(a,+r-)p+ • • • 

X rji  1 j A lip 


1 


1-1  1 

1 * 7" — * -g — ( a i + t — ) p + • • ■ 

A(P) f z ' AB] (p)  1T  fT  1 fT 


lim  1.  I 

p-+0  p * 


A (p) [ z B1 (p)  T 


1 1 Dl 

rr1  = '?r(ai+rr) 


(4-17) 
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liquating  the  right  hand  sides  of  (4-16)  and  (4-17)  results  in 


gT 


N 

r 


L ♦ f-(ai*r) 

i = l f)i^A(-6i)[z  lBK(-p.)  *T  *T 


(4-18) 


This  completes  the  proof  of  liquation  (4-13). 

Q . H . D . 

Lemma  4.3:  Let  the  chain  matrix  parameter  A(s)  of  the  net- 
work of  Figure  4.1  be 


A(s)  = 


Nn 

K.  IT  (s  + z ) 
ai=l  ai 


[da(s)db(s)] 


7T 


and 


A(s) 


(s |na(s)nb(s)=0} 


= 1 


where  -z  t [y (s) z (s) ,k] . Then  the  minimum  value  for  the 
ai 


product  fjgj  is  given  by 


N . 

(fTgT)  - [Z  - r-n-- -]  , (4-19) 

1 1 min  i=l  a. | A't-ou) | 


where  cu  are  the  zeros  of  A(s) 


y(s)z(s)=p. 

The  minimum  of  fjg^  is  achieved  for  the  network  of 
Figure  4.1  whose  short-circuit  output  impedance  is 


Z2,s<«»  = ~ 


z 


1 


y22  /(fTgT)  i = l A^-ap  | (p+aj) 


. (4-20) 


min 


1 


4 
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1 " ‘ 


Consider  < P2  < • * * < 3n  as  adjustable  parameters  of  the 

network,  varying  in  the  intervals  established  in  Equation 
(4-21  ). 

The  problem  is  to  find  for  what  values  of  Bj , P2 , '*'» 

0n  in  the  open  intervals  such  that  the  product  fygy  has  a 

minimum  and  also  to  determine  this  minimum  value. 

First,  observe  that  for  inhomogeneous  ladders  ending 
in  a shunt  arm  admittance,  Bn  * ”,  i.e.,  let 

fTgT  = f(Bj,  f?2,  •••,  6n,1).  (4-24) 

At  the  point  where  the  minimum  occurs, 


Since 


where 


(fTgT)  - * 0 

6K  bk 


(4-25) 


fygy 


1 

aiZA>(-°i)[z"1B1]  (-a4) 


(4-26) 


= f (Bj  , B2»  ’ * ’ * en-l^  * 

i N- 1 

[ z n.  1 (p)  = TT  (1+f  1,  (4-27) 

1 i-1  pi 


expanding  liquation  (4-26)  gives 


rT*T  = 


N 

>. 

i = I 


(•i  A'f-,..) 


I z Bj  | ( *•»  j ) 

(X  . 

(i-/) 

hk 


d-y') 

bk 


wliere  the  bracketed  term  is  independent  of  6^.  Thus 


( fTgT) 


N 

>: 

i = 1 


6 


K 


|2  BjJ (-a.) 

(X  . 

(l-r') 

6k 


ai  2 
d-g-) 

6k 


Simplifying , 


r>tfT*T) 

6K 


N 

E 


dr 


= o 


K i =1  o^A^-a-Hz  B^] (-ai) (-$K+ai) 


(4-28) 


Therefore,  the  numbers  - 6 K ( k = 1,  2,  •••,  n-1)  which 
minimize  f.pgT  = f(Bj,  B2>  *’*»  Bn.j)  satisfy  the  zeros  of 
the  following  rational  function: 


L (p)  = £ 

i = l OjA'f-a^lz  1 Bj ] (-Oj)  (p+Oj) 


(4-29) 


The  poles  of  the  rational  function  L(P)  are  obviously  -a^, 


- a . 


f -a  . Expansion  of 


4 


68 


r 


| z ' 1B1  ] (p)  N [z^BjK-Cj) 

» r 


-l, 


A(pF  “ i=l  JPT^I7TF=IT 


(4-30) 


can  bo  used  to  estimate  l.(p)  within  a constant  multiplier 
Therefore , 


[ z ^B j ] (-c^) 

aiA'(-oi)[z'iB1](-ai)(p+ai)  A' ( a^Hp^T 


“T‘ 


(4-31) 


Since  A'(-a.)|z  ^Bj|(-«j)  > B for  i 
it  follows  from  liquation  (4-30)  that 


= 1 , 2 , • • • , n , 


[ z 1 B j ] (p)  N [z"1B1]  (-ai) 


-1, 


T^Tpl  = ifj  A'  ( -c*i  J (P  + aiT 


N 

Z 


i*l  /Xoi  A' (-«i) (p+«i) 


N 


/X  i=l  /oi|A'(-ai) | (p+a.) 


(4-32) 


Solving  for  x at  p=0  results  in 


[ z ” 1 Bj ] (0)  j N i 

✓X  i-1  ai3/2|A:(-ai)| 


xroy 


i 


* 


69 


Since  [*‘1B11(0)  = 1,  A(0)  = 1,  then 


N , 

J A * £ yyw 

i = l 0i3/Z|A"(-«.) 


(4-33) 


The  short-circuit  output  impedance  can  now  be  written  as 


y22  lPJ  STFI 


t z ” 1B1 ) (P) 


fT  N 

— I 


/A  i *1  /ai  A"(-ot.)  (p+a^) 


(4-34) 


Substituting  Equation  (4-31)  into  Equation  (4-22)  yields 


N 1 

(fTgT)min  = .f  ” T.  T 

l-l  a4  A^(-oi)- 


i-1  a.’ 


|A-(-ai)  | 


'|A'(-a.)| 


This  completes  the  proof  of  Lemma  4.3 


4 


q.h.d. 


Lemma  4.4:  Let  the  chain  matrix  parameter  B(s)  of  the  net- 
work of  Figure  4.2  be 


Ws+V 

z ^sJBCs)  = -2-1 


z 1 (s)B(s) 


[da(s)db(s)] 


[s|na(s)nb(s)  = 0] 


where  -z.  e[y(s)z(s),  k] . Then  the  quantity 


(n+l)/2  . 

8T  = 4 >:  - 

J‘l  822j-l|U'lBl't-62j-1)| 


(4-35) 


where  B ^ are  the  zeros  of  z’  (s)B(s) 


y(s)z(s)  = p. 


gives  the  smallest  possible  value  for  the  total  admittance 


Proof  of  Lemma  4.4:  Let  y(s)z(s)  = p,  then 


l*  BI<P)  fTI(l*JL) 

‘ i-1  “i 


The  parameter  D(p)  = TT  (l+-£-)  is  introduced  with 

i = l ®i 


4 
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6.(i  ■ 1,  2,  n)  arbitrarily  chosen  in  the  intervals 


6 i _ i < 6i  < 6i  i - 1,  2,  •••,  n (4-36) 


Consider  6j,  «2,  •••,  «n  as  parameters  of  the  system.  Use 
Equation  (4-13)  which  is  repeated  here  as  Equation  (4-37), 


gT  " 


N i i bl 

l — . *•  f“(ai  + • (4‘37) 

i=>l  BiZ[z‘1B]'(-8i)D(-Bi)  *T  rT 


This  identity  shows  that  the  total  admittance  gT  can  be 
expressed  as  a rational  function  of  6^,  62,  *’’* 


gT  « F(61,  «2,  • *,  «„)• 


(4-38) 


This  function  has  at  least  one  minimum  inside  the  intervals 
defined  in  (4-36).  At  the  point  where  the  minimum  occurs, 


*K  i*l 


B.[z*1B]"(-Bi)D(-Bi)(-«K+Bi) 


fT6K 


(4-39) 


for  k = 1,  2,  n.  The  derivation  of  (4-39)  is  similar 

to  that  of  Equation  (4-28).  From  Equation  (4-39),  it 
follows  that  the  numbers  -6^  (k  * 1,  2,*’’,  n) , which 

implies  minimizing  gT , satisfy  the  zeros  of 
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L(p) 


N 

£ 

i-1 


6i[z 


B]'('Bi)D(«p.)(p+Bi) 


1 

rip* 


(4-40) 


By  observation  it  follows  that  the  poles  of  L(p)  are 
the  numbers  -Bj,  "6n*  Expansion  of 


yields 


Q(p) 

P [ z ” 1 B ] (p) 


D(p) 

plz^BKp) 

n D C - B±)  l 

i = l Bi[z-:1BJ'(-Bi)(p+Bi)  *1* 


(4-41) 


The  foregoing  expression  should  be  able  to  estimate  L(p) 
within  a constant  multiplier  which  implies 

D2(-Bi)  = 1,  i « 1,  2,  “ ' , n.  (4-42) 


Therefore,  it  follows  from  Equation  (4-41)  that  for  gT  to 
be  a minimum, 

J-  * p t ■ ,T 1 j (4-43) 

[z'1B] (p)  *T  i-1  |(I  B] ' C-B±)| 

If  A (p)  is  used  instead  of  D(p),  an  identical  expression 

is  obtained  for  6-^ Hence,  for  gT  to  be  a minimum, 

[z^BKp) 

A(p)  * D(p) , so  the  ladder  is  symmetrical  as  shown  in 
Figure  4.3.  The  parameters  A* (p)  and  [z_1Bf(p)  of  the 
f i r < t half  of  the  ladder  network  can  be  obtained  as  follows. 


A ( p) 

l z ' 1 B 1 ( p ) 

1 zC] (p) 

D(p) 

A (p) 

| z 1 B | * ( p ) 

IMP) 

1 z * 1 B]  * (p) 

* 

* 

* 

* 

|zC|  (p) 

n (p) 

1 zC 1 (P) 

A (p) 

f z ~ 1 B J (p)  = 2A*  (p)  I z"  1B]  * (p) 


= fT  TT  (1+^-).  (4-44) 

T i=l  Bi 


Since  A (p) 


p=0 


* _ i * 

= 1,  zeros  of  A (p)  and  I z B]  (p)  inter- 


lace in  the  intervals  defined  in  (4-36).  Therefore, 


* (n+l )/2 

A (p)  = TT  (l  + o-0— ), 


j = i 


B2j-1 


(4-45) 


and 


-1  * n 

(z  AB]  Cp)  = ‘ TT  (l+o2-) 

1 J=1  B2j 


(4-46) 


Then  applying  the  identity 


g = TT  — j rr 

i = l 0|A'(-ai)  [z  1B](-ai) 


(4-47) 


for  the  first  half  of  the  ladder,  yields  the  following 


express i on : 
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(n+1 ) / 2 
in  i n i muni  g.^.  = 4 >: 


, — : i • 1 4 - 48 ) 

j-l  s‘2j.,|U  lBI'l-82j.,)l 


This  completes  the  proof  of  l.cmma  4.4 


Q.li.l). 


I:or  the  inhomogeneous  ladder  netwprk  of  Figure  4.1, 
the  open-circuit  voltage  transfer  function  is  given  by 


V2(s>  _ J 
Vj  fs)  ■ A(s) 


(4-49) 


where 


Nn 

Kq  TT  (s+a  ) 


A(s)  = 


ai=l  ai 


(da(s)db(s)] 


TT1 


and 


A(s) 


= 1, 


f s ln_(s)n. (s)  = 0] 


and  -zfl  e (y(s)z(s),  k] . If  y(s)z(s)  = p,  Equation  (4-49) 


becomes 


A (p)  = -jq — — 


(4-50) 


» (p+a,) 

i = l 1 


where 


M = n a . 


i = l 


To  synthesize  the  ladder  network  one  makes  use  of  the 
following  relationship  between  z parameters  and  chain  para- 
meters : 


4 


i 


z,-,(p)  1 /a(p)  , 1 

--^py  - jjcjTpymJTpj  *Tpr 

11 


7 S 


It  is  possible  to  construct  Zjj(p)  as  follows. 


where 


and 


H 


(p)  = 

a(p) 


A(p)  = ka(p) 

[ zCJ (p)  |zCj (pj 

N 

= * (p+o.), 

i = l 1 


N-l 


[ZC] (p)  - P w (P  + Y: ) » 

i*l  1 


(4-51) 


and  0 - Y1  * “l  * Y2  * < Yn  < °n»  (4-52) 

where  y^  are  to  be  determined.  Once  they  are  found  then  the 
continued  fraction  expansion  of  z^^(p)  about  p = °°  yields 
the  inhomogeneous  ladder  network  shown  in  Figure  4.1.  For 
this  purpose  let  the  objective  function  be  defined  as 

F =■  f.p  + WpgT . (4-53) 

Definition  4.5:  Synthesis  of  an  inhomogeneous  ladder  network 
from  a given  transfer  function  is  said  to  be  optimal  if  the 
realization  of  liquation  (4-51)  is  obtained  with  minimum 
value  of  F. 


Lemma  4.5:  I f an  inhomogeneous  ladder  network  is  optimal, 
then  (a)  f^  = w^.g^,  and  (b)  the  solution  for  yj , y2>  Yn 


4 
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is  independent  of  w,, . 

Proof:  The  variable  k in  Equation  (4-51)  is  merely  an 
impedance  scaling  factor.  Therefore,  the  objective  function 
F * fj  * Wpg.p  can  be  written  as 


F = kf„,  + 


T1  + k wFgTl' 


(4-54) 


where  kfT1  * fT  and  ^gT1  = gj.  For  the  network  to  be 
optimal,  the  derivative  of  F must  vanish.  Thus 


3F  - 2 1 

FIc  * fT1  ' k wpgTl  = °* 


Solving  for  k yields 


wFgTi 


(4-55) 


Substituting  the  optimal  value  for  k gives 


fy  = kf^,j 


FgTl 


frp^  =/Wpf,j.jg.p^  9 and 


WpgT  s wp  ^ gT1  gTl  = */WpfT1gT1 . 


Therefore  fT  * WpgT  = ^wFfxi8Tl  Proves  (a) 


To  prove  (b) , note  that  for  the  network  to  be  optimal 
the  derivative  of  F with  respect  to  must  vanish.  Hence 


* -I-  * '■*»’  ■ 


% ■ k T?7  * K"f  -W7’  °>  1 ■ 2>  •"'  "•  (4-56) 


Substituting  liquation  (4-55)  into  liquation  (4-56)  yields 


yyT^fTgT^  = °*  1 " lf  2»  ’ 


(4-57) 


The  solution  of  Equation  (4-57)  for  y^,  y2*  *’’»  Yn  is 
independent  of  both  k and  wp.  This  completes  the  proof. 


Q.E.D. 

Lemma  4.5  makes  it  possible  to  let  wp  = 1 without  any 
loss  of  generality.  Consequently,  it  is  only  necessary  to 
find  the  network  which  minimizes  the  function,  F = fT  + gT, 

Lemma  4.6:  For  the  inhomogeneous  ladder  network  of  Figure 
4.2,  the  minimum  value  of  the  function  F = f^,  + wpg^  is 

obtained  if  and  only  if 


(a)  f,  ■ 


/ i=l ,2 , * • * , jn, 
i fn+2*i  \ i=l ,2 , • • • , l(n+l) ; n 


n even 


isl  2 * * * .n 1 

.11,4,  Jll  , 

«i  = 8n*l-i  t,.!  2 ...  l(n 


n even 


y(n-l) ; n odd 


(b)  fq.  - Wpgj. 


4 
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i 


Proof  of  Necessity: 

(a)  The  open-circuit  admittance  parameter 
network  is  given  by  Equation  (4-51): 


1 

zn(.PJ 


t 


of  the 


y* 


11 


(P) 


1 = tzC] (p) 

z11tp)  TaTpT" 


N 

where  a(p)  = it  (p+o.)-  Thus  an  equivalent  optimization 
i = l 1 

problem  is  to  determine  k and  [zC](p)  so  that  the  realiza- 
tion of  y^jjCp)  yields  the  optimum  network. 

Since  y*22(p)  must  similar  in  form  to  y*11(p),  and 

in  particular,  must  have  the  same  poles,  let 


y*22(p) 


[ zC] * (p) 

J^aTpT 


where  k*  and  [zC]*(p)  are  to  be  determined  so  that  the 
realization  of  y*22^p^  yiel^s  the  optimum  network.  The  two 

optimization  problems  in  terms  of  y^^Cp)  and  y*22(P)  are 

identical  and  must  lead  to  identical  solutions.  Therefore, 


y*n(p)  * y*22(p)  (4-58) 

is  a necessary  condition  for  the  optimum  network.  The 
element  values  are  obtained  by  the  continued  fraction  ex- 
pansion of  y*2j(p)  [or  y*22(p)  or  z j ^ Cp) ] • Since  the 

continued  fraction  expansion  is  a unique  process,  Equation 
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(4-58)  implies  physical  symmetry,  that  is,  condition  (a). 
The  necessity  of  condition  (b)  has  already  been  established 
by  Lemma  4.5. 

Proof  of  Sufficiency:  To  prove  sufficiency,  it  must  be 
shown  that  the  network  which  satisfies  both  conditions  (a) 
and  (b)  and  has  the  prescribed  transfer  function  is  unique. 
This  is  established  by  using  the  fact  that  Figure  4.2  is  a 
compact  network  in  p-domain.  Thus, 


11 


N 

i?i 


(i) 


22 


‘22 


N 

£i  p*“i 


(i) 


where  k^1^  k22^  = [k12^1^]2,  i = 0,  1,  2,  •••,  n and 
kn  = k22  > 0.  This  implies  that  k^  = k22  = |k12|. 

yi,(p)  . . 

can  be  obtained  (to  within  a constant  multiplier)  by 

expanding 


1 

pklz'^Bj  (p) 


in  partial  fractions  and  changing  the  sign  of  any  negative 
residues.  The  continued  fraction  of  the  function  so 
obtained  and  impedance  scaling  to  adjust  fT  * Wpgj  yields 

the  element  values.  The  synthesis  process  is  unique  and  can 
lead  to  only  one  network  for  a given  A(p).  This  completes 
the  proof. 

Q . E . D. 


A 
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From  the  method  described  in  the  proof,  it  follows  that 


1 


1 


PznfpJ  Plz  B](p) 


IT 


P TT  (1+#-) 
i = l Bi 


^o+  " A_ 

p i-i  p+Bi 


where  the  residues 


kn  3 Pi 


Pi  z" 1 B ] (p) 


|p=0 


[ z ~ 1 B ] (0) 


ki 


P [ z ~ 1 B ] '(p) 


N 


1 

"Itii  <•  i 


-B.[z  B] " ( - B i ) 


P=*Bi 


1 . ^ 1 
Tf  + 2_  TT 

An  l f a ^ i _i  r%  r _ li 


pz11(p)  p[z  B] (0)  i = 1 -Bi[z  B] '(-Bi) (p+Bi) 


N 


' 1 m -*  , 


ZU(P)  [z  B]  (0)  i = l -B.[z  lB]  '(-B.)(p+B.) 


1 


(z^BKO) 


N 


1 + 


[ z'2B] (0)p 


i=i  6± | |z"ab] '(-b±) | (p*e±) 

(4-59) 

Lemma  4.7:  If  the  inhomogeneous  ladder  network  of  Figure 
4.1  is  optimal,  then 

f i 3 wp8n+l - i » i 3 1,  2,  •*,  n.  (4-60) 


i 


Proof:  Without  loss  of  generality  let  wp  = 1 , Then 
Equation  (4-60)  becomes 

fi  * 8n+i - i * i = 2»  *•*»  n-  (4-61) 


Equation  (4-61)  exhibits  the  characteristic  of  an  antimetri- 
cal  network  which  can  be  defined  by  the  following  equation. 


zll(p)  = £ y22(p)* 


(4-62) 


For  the  inhomogeneous  ladder  network  of  Figure  4.1, 
it  is  well  known  that 


zn(p) 


fca(p) 

pTzC  j/J  tp)  ’ 


(4-63) 


N-l 

where  [zC1](p)  = E (p+Y-)-  The  continued  fraction  expan - 
1 i = l 1 


sion  of  z11(p)  yields  the  network.  In  order  to  derive 
Equation  (4-62),  an  alternate  approach  using  y parameters 
is  employed.  One  can  write 


X 22(p) 


_MP) 

[z_1B] (p) 


ka(p) 

■ h^bicpj 


ha(p) 

i^Bi-fp) 


(4-64) 


— 
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where 


r z 1 b]  (P) 


N-l 

TT  (l+f-), 

i = 1 pi 


and 


( z" 1B] * (p) 


N-l 

TT  (p+B- ) . 
i = l 1 


The  equivalent  optimization  problem  is  to  determine  h and 
[z'*B]*(p)  so  that  the  realization  of  y^Cp)  yields  the 
optimal  network.  Note  that  the  form  of  y22(P)  *-n  Equation 
(4-64)  differs  from  that  of  Zjj(p)  in  Equation  (4-63)  only 
in  that  y22(p)  has  no  Pole  at  the  origin.  Let 


y*22(p) 


y22(p)  = ha(p) 
p p [ z _1B1 * (p) * 


(4-65) 


The  continued  fraction  expansion  of  y*22(p)  about  p = 00 

* 

yields  the  ladder  network  shown  in  Figure  4.4.  Let  F = 

* * * * 

fj  + gT  * where  f^.  and  gT  are  the  sum  of  impedances  and 

admittances,  respectively,  in  Figure  4.4.  Since  Figure  4.4 

* 

is  the  dual  of  Figure  4.1,  minimization  of  F in  Figure  4.4 
is  equivalent  to  minimizing  F in  Figure  4.1.  Finally,  the 

problem  of  minimizing  F in  terms  of  Zj ^ (p)  given  by  Equation 

* 

(4-63)  is  identical  to  that  of  minimizing  F in  terms  of 
y* 2 2 (P ) J»iven  by  Equation  (4-65).  Therefore,  identical  so- 

- 1 * 

lutions  will  be  obtained,  that  is,  h=k,  and  [z  B]  (p)  = 

[zC  ](p).  Thus  Equation  (4-62)  is  established  for  the 


•j  t-nsr 


r 1 
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«Sl» 


optimal  network  in  the  case  of  Wp  = 1.  This  completes  the 
proof. 

Q.H.D. 

F.xample  4.1:  This  example  illustrates  the  optimal  synthesis 
procedure  of  the  inhomogeneous  ladder  network  that  is  shown 
in  figure  4.1.  Let  the  chain  parameter,  A(s),  of  the  ladder 
be  given  as 


A(s)  = 


s6+23s5+236s4+1338s3+4477s2+8287s+7238 

36(s3+15sZ+7Ss+12S) 


It  is  desired  to  realize  a ladder  with  minimum  total  im- 
pedance-total admittance  product. 

Solution:  The  chain  parameter  A(s)  can  be  factored  to  dis- 
play its  poles  and  zeros, 

(s+2±j  1.732)  (s  + 3. 5±  j 3. 122)  (s+6±j3.317) 

” ' s 1 3 • 


36 (s+5) 


Some  algebraic  manipulations  yield 


A(s) 


^(s  + lHs  + 2) 
s + 5 


^(s  + DCs+2) 
4 (s+5) 


^(s  + 1)  (s  + 2) 
9 (s+5) 


As  can  be  seen  from  Figure  4.5,  the  zeros  of  A(s)  interlace 
with  respect  to 

' (s  + l)(s  + 2)> 


s + 5 


i 


' 
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for  k = 1,4,9. 

The  next  step  in  the  synthesis  process  is  to  pick 
arbitrarily 

*(s)  - and  y (s)  = s+2. 

Hence , 

A(s)  = [l+y(s)z(s)J  [ 1+^-y  (s)  z (s)  ] [ 1+g-y  (s)  z (s)  ] . 


By  use  of  the  mapping  p = y(s)z(s),  A(s)  is  transformed  into 
p-plane  as 

A(p)  - (l+p)(l*|)(l+g). 

The  minimum  value  of  the  total  impedance-total  admittance 
product  is  obtained  from  Equation  (4-19): 


fygy 


where 


> i 

min 

i = l aVZ|A'(-ai)| 

A'(-o,)  = A " ( - 1 ) = t, 


A'(-a2)  = A'  (-4)  = 


A'(-a,)  = A ' (-  9)  = 


and 

f-pgj 


min 


372TTT 


12 


+ T7T.'*"-7.  + T7TT*T 


<*1  |A'(-ax)|  oj'  IA'C-OjJI  otj  4 1 A ' (-a^)  | 


(U)2 

t“5_)  ' 


4 


Therefore,  for  optimal  ladder,  fT  = = -g-. 

This  minimum  value  is  achieved  by  the  ladder  with  the 
short-circuit  output  impedance  given  by  liquation  (4-20). 

_i__  = [r l ... .j 

y22^  /(fTBT)  i = l (p  + ou) 

min 


j(p+l)  2(^-)(P+4)  3(^)(p+9) 

|(p2*13p+36)  + |(p2+10p+9)  + ^(p2  + 5p+4) 
(p+1) (p+4) (p+9) 


Thus , 


3p2_+33p+6(> 

+ 1 4p^+49p+36 


22 


(P) 


p5+14p2+49p+36 

3p2+33p+66 


The  optimum  values  of  impedance  level  and  admittance  level 
arc  shown  in  Table  4.1.  The  corresponding  ladder  network 
is  shown  in  Figure  4.6. 

lixample  4.2:  This  example  demonstrates  the  synthesis  pro- 
cedure of  the  optimal  inhomogeneous  ladder  network  of 
Figure  4.2  with  minimum  total  admittance  level,  gy.  Let 


the  chain  parameter,  B(s),  of  the  network  be  given  as 


AD-A064  291 


UNCLASSIFIED 


NAVAL  ACADEMY  ANNAPOLIS  MD  DIV  OF  ENGINEERING  AND  WEAPONS  F/G  12/2 
SYNTHESIS  OF  OPTIMAL  LADDER  NETWORKS. (U) 

NOV  78  T S LIM 

USNA-EW-16-78  NL 


Table  4.1  Optimal  Network  of  Example  4.1 


f3  = 1 


Bt  * 1 


1 

14 

49  36 

3 

33 

66 

3 

27 

36 

6 

30 

12 

36 

12 

36 

nr  o 4(s7+19s6+158s5+738s4+2074s3+3470s2+3149s+1155) 

t -7 

3(s°+12s^+48s+64) 


Solution:  The  chain  parameter  B(s)  can  be  factored  to  dis- 

« 

play  its  poles  and  zeros, 


B(s)  = »(s  + l)  (s  + 2.5*j0.866)  (s  + 3±jl  .414)  (s  + 3 . 5± j 1 . 658) 

6CS  + 4)"5 


Some  algebraic  manipulations  yield 
B(s)  = 8(s*l)  fl 


_ (s  + 1)  (s+3)1  L (s  + 1)  (s  + 3)l  L(s*l)  (s  + 3) 

1 — s+T 7T**Tj — ~3Ts+TJ — • 


4 


*■  * >*■%■& 
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Arbitrarily  picking  z(s)  = s + 1 and  y(s)  - 1^,  results  in 


z'^sJBts)  - 8[l+y(s)z(s)]  [1  ♦yy (s) z (s)  ] [ 1+j-y  (s)  z (s)  ] . 


As  can  be  seen  from  Figure  4.7,  the  zeros  of  z_1(s)B(s) 
interlace  with  respect  to 


(s  + lHs  + S),  k 
s+4 

W J 


for  k = 1,2,3.  The  mapping  p = y(s)z(s)  yields 


[z_1B]  (p)  » 8(l+p)(H-£)(l+£). 

The  minimum  total  admittance  level,  gT,  is  obtained  from 
Equation  (4-35): 


tz  B]  C-0j)  33£|[z‘aB]'(-B3)| 


where  Bj  - 1 , B2  = 2,  B3  = 3,  [z"1B]  '(-Bj)  = f, 
[z  B)  (-e3)  -j,  and  g?  = 4[g^y  + 
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The  ladder  with  the  specified  [z_1B](p)  that  achieves 
this  minimum  total  admittance  level  is  a symmetrical  ladder 

A _ l A 

as  shown  in  Figure  4.3.  Let  A (p)  and  [z  B]  (p)  be  the 
corresponding  parameters  of  the  first  half  of  the  symmetri- 
cal ladder.  Then  it  follows  from  Equation  (4-44)  that 

[z  *B]  (p)  = 2A* (p)  [z_1B] * (p) 


- 8(l+p)(l+f)(l+J). 


The  short-circuit  output  impedance  — of  the  first  half 

y (P) 

of  the  two-port  network  in  Figure  4.3  is 


1 

y*(p) 


Z*(P) 


u-ytp)  _ 4(1^ 

A (p)  (l+p)(l+|) 


Thus 


y*(p) 


p^+4p+3 

6p+12 


The  values  of  impedance  level  and  admittance  level  are 
listed  in  Table  4.2.  The  corresponding  ladder  network  is 
shown  in  Figure  4.8. 
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tr- 


iable 4.2  Half  of  Symmetrical  Ladder  of  Example  4.2 


1 

S 


1 4 3 

6 12 


f 


2 


3 


*1 


2 

I 


2 3 

3 


f 


1 


= 1 


3 


Table  4 . 3 Element  Values  of  Example  4.1 


**•1 

0.333 

! Rzai 

< 

0.500 

RT 

3aj 

1.000 

Rla2 

0.067 

' R2a 

j 

0.100 

1 

| 

R3a2 

0.200 

Lla 

0.067 

1 L2a 

j 

0.100 

L3a 

0.200 

Cla 

3.000 

! C2a 

i 

2.000  j 

C3a 

| 

1.000 

Rlb 

i 

0.500 

— 

R2b 

1 

.000 

1 

! R3b 

1.500 

Clb 

1.000 

C2b 

0 

.500 

C3b 

0.333 

gT 

= 

fT 

= 

11 

6 

min 

min 
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Table  4.4  Element  Values  of  Example  4.2 


= i 
3 

min 


[ 
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CHAPTER  S 


TIME  DOMAIN  RESPONSES  OF  A CLASS  OF  OPTIMAL 
RC  LADDER  NETWORKS 

Section  5 . 1 Introduction 

This  chapter  is  concerned  with  the  time  domain  responses 
of  certain  classes  of  optimal  ladder  networks  as  shown  in 
Figure  5.1.  The  state  equation  of  the  network  is 
[V]  = [A][V]+[B]u  where  the  state  variable  [V]  represents 
the  voltages  across  the  capacitors.  It  was  found  that  if 
the  eigenvalues  of  the  tridiagonal  matrix  [A]  are  prescribed 
in  a certain  pattern,  then  the  corresponding  optimal  ladder 
can  be  obtained  systeinmatically  by  using  the  formulas  pre- 
sented in  the  theorems  that  follow. 

Section  5.2  Synt hesis  Procedures 
Definition  5.1  Let  a tridiagonal  matrix  be 


If  the  off-diagonal  elements  of  [A]  satisfy  the  condition 
that  a^  = a(n-i  + l) (n- j + 1)  ’ then  thc  tridiaRonal  matrix  [A] 
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is  said  to  be  antimetrical . 

Consider  the  doub le-terminated  RC  ladder  network  of 
Figure  5.1.  Let  the  state  equation  be  fV]  = [Al(V]+[Blu 
where  the  state  variable  [V]  represents  the  voltages  across 
the  capacitors.  The  |A)  matrix  is  tridiagonal  as  shown 
below . 


1 f 1 + 1 l 
L1  K1  r2 


R2CZ 


- 1 r 1 + 1 1 


R3^\3 


C3  U3  R4 


rt 

n n - 1 


1 - 1 r 1 + 1 i 

nr  c-  r~7 

n n n n n+1 


(5-1) 


Let  an  RC  ladder  network  be  double- terminated  as  shown  in 


Figure  5.1  where  Rj  = Rn+2* 
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n(n-l) • • • (n-k+2)Rj 
(k'l) I 


(5-2) 


where  k * 2,3,* • • ,n+l , and. 


2(k-l)l 


n(n-l) • • • (n-k*l)p1R1 


(5-3) 


where  k = l,2,***,n,  [Observation:  Equations  (5-2)  and 
(5-3)  imply  that  Rk  = Rn.k+2  where  k * 2,3,***fn+l,  and 

Ck  * Cn-k*l  w*iere  ^ = l,2,**-,n.]  then  the  impulse  response 


is 


N -P;t 

V‘>  - kic 

i=l 


where  p^  * ip^ 
and 


kj  * (-1) 


i + 1 


H 


(i-1)! (n-i)!Pln_1 


(5-4) 


(5-5) 


and 


R . . N 

n + 1 -it 


i = l 


N 

TT  R 7TC. 
i*l  1 i = l 1 


(5-6) 


For  the  n=2  case,  the  voltage  transfer  function  is 
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v») 

Vsfs) 


2W 

~2 — 1 2 

s *Irr(:rs*H12ciz 


rici 


Rici 


2Ri2c7  s+j quj  s^7 


(5-7) 


Lquation  (5-7)  shows  that  the  impulse  response  is 


VQ(t)  = kje  Pjt  + k2e  P£t 


where 


P1  = R'VjT'  p2  = R^,  = 2*V 


kI  = RlCiM»  k2  = RiciH» 


2 — 7' 

2Rj  Cj 


lor  the  n=3  case,  the  voltage  transfer  function 


V.!J 

vs(7) 


4 . SRj  Cj'  s5^18R]~3C1^s2  + 22R1^C1s  + 8R1  * 


9R  ^ 
y 1 L1 

r *J  * - 4Z4  • 

Rll'l  SRj^Cj2  9Rj!C^ 
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_ 

9 R j 3Cl3 


( s + JTT[T^)  ( s + SU'JC x}  ( s 4 K|Tj } 


3 — J 

9RX  CL 


2r  2 2r  2 ^r  2r  2 

ffRl  L1  TK1  L1  . ffRl  C1 

2 3 — 2 

s+3T^uj  s+nqrY  s+iq^ 


(5-8) 


Equation  (5-8)  shows  that  the  impulse  response  is 


-Pit  -p2t  -p,t 

VQ(t)  = kje  1 +k2e  L +k3e  J , 


where 


P1  = W^> 


p2  = ZW^  = 2pl’ 


p3  ~ r^y  = 3pi 


and 


k - ^R  2r  2h 
ki  fri  li  h» 


Ir^Cj2", 


k,  = 


9U  2(,  2 

8'l  l'l  n» 


and 


H = 


3 — S’ 

9Rj  Cj 


Values  of  these  two  cases  agree  with  the  results 


specified  in  liquations  (5-4),  (5-5),  and  (5-0). 

It  is  noted  that  the  value  of  H given  in  (5-0)  is 
exactly  the  maximum  gain  of  Kuh's  optimal  ladder  (12j. 

Lemma  5.1:  Let  an  KC  ladder  network  he  doubl c- term i nated 

as  shown  in  figure  5.1  where  Rj  = l<n+j.  If 


n(n-l) • • • (n-k+2)Rj 
(k-l)l 


(5-9) 


where  k = 2,3,***,n+l,  and 


2 (k - 1 ) ! ^ 

n ( n - 1 ) * • • ( n - k + 1 ) p j 11  j 


(5-10) 


where  k = L,2,***,n,  then 

(a)  The  diagonal  elements  of  the  [A]  matrix  of  (5-1) 
arc  equal  and  each  is  equal  to  the  sum  of  the  eigenvalues 
divided  by  the  number  of  eigenvalues;  i.e.. 


1 - 1 ) 


(n+l)Pl 

2 


(5-11) 


(b)  The  | A]  matrix  is  antimetrical  as  defined  in 
lie  fi  n i t i on  5.1. 

Proo  f : 

(a)  Since  pj  = ipj  for  n = 2,3,*'',n,  therefore,  the 
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n 

^Pi 
i = l 1 


n(n+l) 

2 P1 


(5-12) 


It  follows  from  liquation  (S-il)  that 


1 ♦ 1 


Rk  Rk*l 


(k- 1)  ! 


k! 


n"(n-l)(n-2').'..  (n-k+ITRJ  + nTn-1)  (n-2)  • (n-k+l)!^ 


(k-1) ! (n-k+l)  + (k-l) Ik 
n(n-l)  (n-2)  • • • (n-k+D  (n-k+l)R~ 


(k- 1) ! (n+1) 

n^n-1)  (n-2)  • ^7n-T<  + 2)  (n-k+l)R 


From  liquation  (5-10)  and  (5-13),  it  follows  that 


1 ( 1 + 1 , 
Lk  Kk  Kk*l 


n(n-l) (n-2) • • • (n-k  + l)pj  R1 

rnc--i)i 


(5-13) 


(k-1) ! (n+1) 

n(n- i ) (n- Z) • • •Tn-V+2) (n-k+TTRJ 


This  results  in  the  equation 


1 f 1 + 1 , 
Lk  Kk  Rk-1 


n+1 
2 P1 


4 
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and  therefore  by  Equation  (5-12),  this  value  is  the  sum  of 
divided  by  n.  This  completes  the  proof  of  (a). 

To  prove  that  the  [A]  matrix  is  antimetrical , it 
suffices,  according  to  Definition  5.1  and  Equation  (5-1), 
to  show  that 


RkCk  Rn-k+2Cn-k+l 


(5-14) 


It  follows  from  Equation  (5-9)  that 


n-k  + 2 


n(n-l) • •• (n-k+2) (n-k+1) • • • kR^ 
(n-k+1) (n-k) •• *k* (k-1) ! 


n(n-l) * • • (n-k+2)Rj 
(k-1)! 


Similarly,  it  can  be  shown  that  CR .^+j  = C^.  Therefore, 


(5-14)  follows. 


Q.H.D. 


Lemma  5.2:  Let  an  optimal  RC  ladder  network  with  minimum 
total  capacitance  and  with  Rj  * Rn+1  be  double- terminated  as 

shown  in  Figure  5.1.  Let  its  impulse  response  be 


n -p.t 
V It)  = E k.e 


4 
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where 


Pi  - iPr 


(5-15) 


and 


kj  “ (-1) 


i + 1 


H 


(i-1) ! (n-i) !pj 


?TT 


(5-16) 


H = 


Vi  " 
55T-  J.Pi 

5“  R 11 
i-1 1 


n n 
TT  R.  TT  C. 
i-1  1 i-1  1 


(5-17) 


Then 


n(n-l) • • • (n-k+2) R, 

R = i, 

(k- 1) ! 


(5-18) 


where  k = 2, 3, • • • ,n+l , and 


2(k-l)! 


n(n-l) • • • (n-k+l)p1R1 


(5-19) 


where  k = 1,2,* • • ,n. 

From  the  given  impulse  response,  it  follows  that  the 
voltage  transfer  function  of  the  ladder  network  is 


Vs) 

Vs(s) 


H 


n 


TT  (s+p.) 
i-1  1 


For  n=3,  the  optimal  RC  ladder  has  the  short  circuit  input 
impedance 
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The  residues  as  expressed  in  Equation  (5-16)  are 

tabulated  in  Table  5.2.  A common  factor  is  taken  out  of  the 
residues  to  make  more  apparent  their  variations.  Plots  of 
impulse  responses  are  shown  in  Figures  5.3  and  5.10. 

Table  5.3  shows  the  residues  of  the  unit  step  responses 
of  optimal  ladders  with  eigenvalues  pi  = ip^.  Figure  5.11 

shows  plots  of  the  responses  of  these  ladders  to  unit  steps. 
Therefore,  a circuit  whose  desired  impulse  response  or  step 
response  has  similar  variations  can  be  realized  as  an  opti- 
mal ladder  by  the  use  of  Equations  (5-2)  and  (5-3). 

Table  5.2  Residues  of  Impulse  Response, 

n i+i  'Pi1 

V_(t)  = Z (-D  . 

i = l * 

of  Optimal  Ladders  where  p^  = ipj 


n 

Common 

Factor 

ki 

k2 

k3 

k4 

k5 

k6 

** 

^1 

Xs 

00 

2 

7P1 

1 

-1 

3 

hi 

1 

-2 

1 

4 

1 

-3 

3 

-1 

5 

^pi 

1 

-4 

6 

-4 

1 

6 

T7?l 

1 

-5 

10 

-10 

5 

-1 

7 

7 „ 

mpi 

1 

-6 

15 

-20 

15 

-6 

1 

8 

TZPl 

1 

-7 

21 

-35 

35 

-21 

7 -1 
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Table  5.3  Residues  of  Unit  Step  Response, 

n i 'Pi1 
V_(t)  = kn  ♦ Z C-l)  k.e 
° 0 i = l 1 

of  Optimal  Ladders  where  = ipj 


n 

Common 

Factor 

k 

o 

kl 

k2 

k3 

k4  k5  k6 

2 

TP1 

1 

-2 

1 

3 

1P1 

1 

-3 

3 

-1 

4 

T6p1 

1 

-4 

6 

-4 

1 

5 

12  pl 

1 

-5 

10 

-10 

5 -1 

6 

6TP1 

1 

-6 

15 

-20 

15  -6  1 

Theorems  5.1  and  5.2  are  concerned  with  optimal  ladder 
networks  where  the  eigenvalues  of  the  [A]  matrix  are  p^=ip^. 

2 

If  instead  the  eigenvalues  are  p^i  pj  with  ki  given  in 

Table  5.4  and  Table  5.5  for  impulse  and  step  response 
respectively,  then  the  element  values  of  the  optimal  network 
can  be  calculated  from  the  following  two  formulas. 


R . 1-3-5-  ••  (2k-3) 

k (2n-l) (2n-3) • • • (2n-2k+3) 

where  k = 2 ,3 , • • • ,n  + l and 


n (n  - 1 ) • • • (n -k  + 2) Rj 
(k- 1)  ! 

(5-23) 


4 
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c = (2n-l)  (2n-3)  • • - (2n-2k+3)  . 2(k-l)  ! 

k 1 • 3 • 5 • ■ • (2k-l)  n(n-l) • • • (n-k+l)p1R1> 

(5-24) 

where  k = 1 ,2 , • • • ,n. 

The  impulse  responses  of  optimal  ladders  with  eigen- 
2 

values  pi  = i pj  are  shown  in  Figure  5.12.  The  responses 
to  unit  steps  are  shown  in  Figure  5.13. 

Table  5.4  Residues  of  Impulse  Response, 

VQ(t)  = £ (-l)i-1kie  Pi\ 

_ i = l 

2 

of  Optimal  Ladders  where  pi  = i p^ 


n 

Common 

Factor 

kl 

k2 

k3 

k4 

k5  k6 

2 

7P1 

1 

-1 

3 

37pl 

5 

-8 

3 

4 

lTpl 

7 

-14 

9 

-2 

5 

5T2?1 

42 

-96 

81 

-32 

5 

6 

5T7pl 

66 

-165 

165 

-88 

25  -3 

The  formulas  derived  and  presented  so  far  in  this 
chapter  provide  straightforward  methods  of  synthesizing 
optimal  ladder  networks  whose  eigenvalues  follow  certain 
patterns.  When  the  eigenvalues  are  not  arranged  in  these 


4 


b 


< 
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Table  5.5  Residues  of  Unit  Step  Response, 

n i -p4t 

V0(t)  - k * 2K-l)Ve  1 , 
i = l 


2 

of  Optimal  Ladders  where  p^  = i 


n 

Common 

Factor 

k 

0 

ki 

k2 

k3 

k4  k5 

2 

SP1 

3 

-4 

1 

3 

T^l 

10 

-15 

6 

-1 

4 

mpi 

35 

-56 

28 

-8 

1 

5 

5T7pl 

126 

-210 

120 

-45 

10  -1 

patterns,  the  existing  synthesis  methods  require  a continued 
fraction  expansion.  The  complexity  of  such  a procedure 
increases  rapidly  with  the  number  of  ladder  sections. 
However,  for  a two-section  ladder  it  is  possible  to  derive 
simple  and  useful  formulas  based  on  the  following  equations: 


1 ( 1 + l , pl+p2 
C7(r7  2 ’ 


C 5 - 2 5 ) 


and 


plp2 


R1*R2*R3 

R1R2R3C1C2 


(5-26) 


Since  Rj  » 
(5-26)  for 


R3  * ■ C2,  solving  Equations  (5-25) 

&2  an<^  yields  the  following  two  formulas 


and 


r' 
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2p. 

R,  “ (5-27) 

P2-Pi 


C.  = — . (5-28) 

Pi 


These  formulas  are  easy  to  apply  and  are  useful  for  two- 
section  optimal  ladder  with  any  two  eigenvalues.  Example 
5.3  illustrates  the  process. 

For  the  n=3  case,  the  short  circuit  input  impedance  of 
the  optimal  ladder  is 

3 2 

. s +a-s  +a.s+a„ 

1 2 1 o 

yll  s3+b7sZ+b.s+b  ’ 

L 1 O 

where 

ao  = P1P2P3 ' 

al  = plp2+plp3+p2p3» 

a2  = pl+p2+p3  * 

u _ P1P2P3P21P31P32 

bo  * 

plp2p21+plp3p31+p2p3p32+p21p31p32 

b = pl2p22p21^pl2p32p31+p22p32p32^p21p31p32(plp2+plp3+p2p3) 
P1P2P21+P1P3P31+P2P3P32+P21P31P32 

L 


* 


Ill 


b = pip2P21(Pl^2^PlP3P31(Pl^3^P2P3p32(P2>P3) 
2 PlP2P21+PlP3P31+P2P3P32+P21p3lP32 

+ P2lP3lP32(Pl't'P2'l'p3j  _ 

P1P2P21+P1P3P31+P2P3P32+P21P31P32 

P21  = P2‘P1 » 

= P3 " pi  > snd 


Example  5.1:  This  example  illustrates  the  optimal  realiza- 
tion technique  for  the  ladder  network  of  Figure  5.1  with  the 


4 
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unit  impulse  response, 


VQ(t)  “ ^r(e't-4e*2t*6e'3t-4e'4t+e'5t) 


The  terminating  resistors  are  Rj  ■ Rfi  » in.  For  n*5, 
Equations  (5-2)  and  (S-3)  yield  the  following  element 


values : 


R2  = Rs  = 5Rj  ■ 5n, 


Rj  ■ R4  - 10Rj  ■>  ion. 


C1  = c5  " Spfrq-  " lF» 


c2  * c4  " iTpjRj  3 nrF>  and 

C3  = = irF- 

The  optimal  network  is  shown  in  Figure  5.14.  The  maximum 
gain  of  the  optimal  ladder  network  is  given  by  Equation  (5-6) 


R,  5 

-T- 

Z R,  11 

i“l  1 


R6P1P2P3P4P5 

R1+R2+R3+WR6 


4 
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15 

= T 

The  optimal  network  is  shown  in  Figure  5.14. 

Example  5.2:  This  example  demonstrates  the  optimal  synthesis 

procedure  for  the  ladder  network  of  Figure  5.1  with  eigen- 
2 

values  p.  = i pj.  The  unit  impulse  response  is 
Vo(t)  = I77(5e't-8e'4t^-3e‘9t). 

The  terminating  resistors  are  R^  = R^  = lfl.  The  problem  is 
to  realize  the  network  with  maximum  possible  gain  H. 

Solution:  Application  of  liquations  (5-23)  and  (5-24)  with 
n=3  yields  the  following  values. 

^2  = R3  = 1>^  * 

Cj  3 Cj  = jF,  and 


The  maximum  gain  H is  obtained  from  Equation  (5-6), 

tI  R4  P1P2P3  45 
M = 

Rl+R2+R3+R4 


The  optimal  network  is  shown  in  Figure  5.15. 
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Example  5.3:  This  example  illustrates  the  optimal  synthesi; 
procedure  for  the  ladder  network  of  Figure  5.1.  The  unit 
impulse  response  is 

Vo(t)  * H(ye" 2t  - ^'5t). 


The  problem  is  to  realize  a ladder  network  with 
possible  H. 


maximum 


Solution:  Substitution  of  the  eigenvalues  into  formulas 
(5-27)  and  (5-28)  yields  the  following: 


t « « 4 

2 7T- 

p2'pl 


c = J_  - 1 
C1  PX  " T 


The  value  of  H is  obtained  from  Equation  (5-6) 


Rl+R2+R3 


iP?  = 3 


The  optimal  network  is  shown  in  Figure  5.16. 
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CHAPTER  6 

CONCLUSIONS  AND  RECOMMENDATIONS  FOR' FUTURE  RESEARCH 


In  Chapter  2 and  3,  theorems  have  been  proved  and 
methods  have  been  demonstrated  to  synthesize  various  types 
of  filters  by  chain  matrix  decomposition  techniques. 

Lossless  (LC)  ladder  networks  terminated  in  a load  resistance 
as  well  as  those  resistively  terminated  at  both  ends  have 
been  investigated.  By  the  use  of  simple  algorithms  pre- 
sented in  Section  2.4,  implemented  by  digital  computer  pro- 
grams, filter  design  can  be  accomplished  with  little 
expenditure  of  time.  Classical  filters  such  as  Butterworth, 
Chebyshev,  and  Bessel  have  been  realized  by  means  of  chain 
matrix  decomposition  techniques.  The  results  of  these 
filters  have  been  compiled  into  tables  for  both  single  and 
double- terminated  cases. 

Chapter  4 deals  with  optimal  synthesis  of  inhomogeneous 
ladder  networks.  It  has  been  found  that  an  optimal 
inhomogeneous  ladder  ending  in  a series  arm  is  symmetrical 
and  the  one  ending  in  a shunt  arm  is  antimetrical . Pro- 
cedures have  been  derived  to  synthesize  these  ladders. 

Chapter  5 is  concerned  with  time  domain  responses  of 
a class  of  optimal  RC  ladder  networks.  The  outcome  provides 
an  insight  into  some  very  useful  properties  of  the  tri- 
diagonal  matrix  [A]  of  the  state  equation  [VJ  = [A][V]+[B]u, 
which  describes  a certain  type  of  double-terminated  ladder 
networks.  It  has  been  found  that  the  [A]  matrix  of  the 
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optimal  ladder  is  antimetrical  as  established  in  Definition 
5.1  and  the  sum  of  the  diagonal  elements  is  equal  to  the  sum 
of  the  eigenvalues  of  [A].  If,  in  particular,  the 
eigenvalues  p^  are  distributed  in  such  a way  that  p^  * ipj  , 

then  the  diagonal  elements  of  [A]  are  equal  and  each  is 

equal  to  the  sum  of  eigenvalues  of  [A1  divided  by  the  number 

of  eigenvalues.  Several  formulas  have  been  found  to 

facilitate  the  calculation  of  element  values  of  optimal 

networks  of  this  type  as  well  as  those  whose  eigenvalues 

2 

have  the  property  of  p-  = i Pj.  Realization  of  these  types 

of  optimal  ladders  is  reduced  to  mere  substitution  of  the 
given  eigenvalues  into  these  formulas. 

While  the  present  thesis  has  been  limited  to  optimal 
synthesis  of  networks  whose  eigenvalues  have  certain  special 
types  of  distribution  pattern,  future  research  should 
include  problems  whose  eigenvalues  have  arbitrary  locations. 
This  thesis  has  investigated  this  type  of  ladder  for  the 
n=2  case  and  formulas  have  been  found  that  greatly  simplify 
the  realization  process.  It  is  believed  that  perhaps 
formulas  can  be  found  for  the  n=3  case,  and  for  the  n=4 
case,  etc.  Sensitivity  analysis  of  these  optimal  ladders 
is  also  an  interesting  topic  for  future  research. 

A worthwhile  topic  in  connection  with  chain  matrix 
decomposition  is  the  derivation  of  the  necessary  and 
sufficient  conditions  for  the  synthesis  of  general  ladder 
network  where  Zj(s)  f z^ (s)  and  y^(s)  f y^ (s) . 

The  ground  work  of  time  domain  analysis  was  laid  in 


4 
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Chapter  5.  Some  distinct  properties  of  the  [A]  matrix 
were  found;  however,  no  more  synthesis  methods  were  found  to 
be  of  value.  Further  study  may  lead  to  a promising  approach 
to  the  time  domain  synthesis. 
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